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I. 


INTRODUCTION 


The  primary  objective  of  this  program  is  to  devise,  through 
the  research  to  be  discussed,  an  optimum  set  of  techniques  ca- 
pable of  compensating  the  residual  boron  always  found  in  detector 
grade  silicon  by  the  process  of  nuclear  transmutation  doping. 

The  material  to  be  developed  by  this  program  is  intended  for  use 
in  infrared  detector  systems  such  as  those  of  the  Air  Force  LADIR 
(Low  Cost  Arrays  for  the  Detection  of  Infrared)  program. 

The  best  silicon  available  today  contains  residual  boron,  a 

12  3 

shallow  acceptor,  in  concentrations  of  the  order  of  10  B/cm  . 

This  is  because  the  segregation  coefficient  of  this  impurity  is 
close  to  unity,  and,  therefore,  will  not  be  substantially  re- 
moved by  successive  passes  in  a zone  refiner.  This  residual 
boron  concentration  has  several  adverse  effects  on  the  operation 
of  extrinsic  photoconductor  IR  detectors  fabricated  from  Ga  or 
In  doped  silicon.  Since  the  B ionization  energy  is  less  than 
those  of  Ga  or  In,  the  B levels  thermally  ionize  at  temperatures 
below  those  of  the  impurities  of  interest  forcing  the  detectivity 
out  of  the  background  limited  condition  at  prematurely  low 
temperatures.  This  effect  requires  additional  detector  array 
cooling,  which  results  in  increased  system  weight  and  cost,  to 
hold  the  detectors  in  the  background  limited  region.  As  a 
result  system  size  and  cost  increase  unless  these  boron 
levels  can  be  populated  by  carriers  from  compensating  donor 
levels.  Furthermore,  the  shallow  boron  optical  transition 
is  beyond  the  atmospheric  transmission  windows  through  which 
the  IR  data  information  is  obtained.  Therefore,  these  boron 

L 
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optical  transitions  can  respond  to  undesirable  background  thermal 
radiation  to  produce  a type  of  nonuseful  background  noise. 

Both  of  these  problems  can  be  eliminated  with  a minimum 
degradation  in  theoretical  detector  gain  if  nearly  exact  com- 
pensation of  the  boron  acceptors  is  accomplished  by  the  addition 
of  shallow  donors.  Greatly  overcompensating  the  boron,  however, 
is  undesirable  because  the  additional  donor  carriers  fill  the  Ga 
or  In  levels  once  all  the  B levels  are  compensated.  This  reduces 
detector  gain  by  reducing  the  number  of  desirable  optical  tran- 
sitions from  which  signal  information  is  obtained.  Additionally, 
the  extra  impurities  from  over-compensation  represent  additional 
charged  impurity  scattering  centers  reducing  the  carrier  mobility, 
which  is  dominated  by  charged  impurity  scattering  at  the  low 
temperatures  at  which  these  detectors  operate.  This  lowering  of 
carrier  mobility  reduces  the  gain  of  the  detector,  which  is  pro- 
portional to  the  product  of  mobility  times  lifetime.  Carrier 
lifetime  is  also  reduced  by  overcompensation  since  additional 
recombination  centers  are  created. 

Nuclear  transmutation  doping  of  silicon  appears  to  offer 
several  advantages  over  conventional  doping  techniques  for  pro- 
ducing the  desired  compensating  donor  concentration.  Briefly, 

the  transmutation  doping  process  is  based  on  the  conversion  of 
30  31 

the  isotope  Si  to  Si  through  an  (n,  y)  reaction  induced  by 

thermal  neutron  capture  during  reactor  irradiation.  The  ^Si 

31 

isotope  is  unstable  and  decays  through  6 emission  to  P,  the 
desired  donor  impurity.  Since  the  production  rate  of  is 
3.355  ppb/10  n/cm  , the  production  of  donors  is  slow  enough  to 
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provide  careful  control  of  the  dopant  concentration.  Further- 
more, the  donors  added  by  transmutation  doping  are  distributed 
uniformly  throughout  the  lattice  because  the  microscopic  abun- 
dance ratios  of  the  Si  isotopes  are  uniform  and  the  thermal 
neutron  diffusion  length  is  very  long  compared  to  typical  ingot 
dimensions . 

The  two  major  advantages  of  NTD-Si,  tight  control  over  doping 
concentration  and  uniformity  of  dopant  distribution,  are  some- 
what offset  by  the  radiation  damage  effects  which  occur  during 
irradiation.  These  effects  cause  no  difficulty  for  float  zone 
Si  neutron  doped  to  a resistivity  of  100  ft-cm  or  lower.  For 
higher  final  resistivity  material,  the  success  or  failure  depends 
intimately  on  residual  impurity  defect  interactions  and  their 
annealing  properties.  One  of  the  major  goals  of  the  present 
research  is  to  define  the  limits  of  final  resistivity  which  can 
be  achieved  relative  to  the  properties  of  the  starting  material. 

The  balance  of  this  report  is  contained  in  five  major 
sections  as  follows: 

• Section  II  develops  the  theoretical  tools  to  understand  the 
expected  electrical  properties  resulting  from  transmutation 
doped  Si  as  functions  of  neutron  fluence.  A discussion  of 
the  transmutation  doping  process  is  also  presented. 

• Section  III  is  a description  of  the  experimental  techniques 
and  apparatus  employed  throughout  this  work. 

• Section  IV  gives  a description  of  the  experimental  results 
and  a discussion  of  possible  interpretations. 
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• Section  V discusses  the  significant  accomplishments 

achieved  during  the  program  and  discusses  areas  of  difficulty 
where  further  research  will  be  required. 

II.  THEORETICAL  ANALYSIS 

The  basic  physics  of  transmutation  doping  in  Si  has  been 
discussed  by  Tanenbaum  and  Mills.1  Although  they  discuss  all 
possible  nuclear  reactions  with  all  the  silicon  and  oxygen 
isotopes  such  as  (n,  y)  , (n,  2n) , (n,  p)  , and  (n,  a) , we  will 
confine  our  attention  to  the  (n,  y)  reactions  with  Si  as  being 
the  only  significant  reactions  in  a thermal  neutron  flux.  This 
viewpoint  has  been  verified  experimentally.1 

The  natural  abundance,  thermal  neutron  (n,  y)  absorption 
cross  section  and  half-lives  of  the  relevant  isotopes  of 
silicon  and  phosphorus  are  listed  in  Table  1. 


Table  1.  Nuclear  Data. 


Natural 

Abundance 

Isotope 

o(n,  y) 

Tl/2 

92.21% 

28Si 

0.08  + 0.03  b 

stable 

b.68% 

29Si 

0.28  + 0.09  b 

stable 

3.05% 

3°Si 

0.11  + 0.01  b 

2.62h 

100% 

31P 

0.20  + 0.02  b 

lU.3d 

At.  Wgt.  Si  = 

28.086,  N(Si 

22  3 

= 4.996  x 10  atoms/cm 

Density  Si  = 

2.33  gm/cm3 

The  nuclear  reactions  of  interest  are  the  following: 


28Si 

(n. 

y) 

29c  . 

Si 

(1) 

29Si 

(n, 

y) 

30c  . 

Si 

(2) 

4 


(3) 


30Si  (n,  y)  31Si  -*■  31P  + B' 


T1/2  = 2 . 62h 


31  32  32 

J P (n,  y)  JZP  - ^S  + B 


(4) 


I1/2  = 14.3d 


The  first  two  reactions  redistribute  the  natural  abundance 

of  the  silicon  isotopes  but  do  not  lead  to  any  doping  effects  or 

radioactivity.  Reaction  (2)  is  quite  fortunate  in  that  it  serves 
30 

to  replace  the  Si  isotope  approximately  as  it  is  burned  up  in 

reaction  (3)  which  produces  the  dopant  3^P.  Reaction  (4)  trans- 

31  32 

mutes  the  stable  isotope  P to  the  unstable  isotope  P which 
32 

B decays  to  S.  Since  reaction  (4)  is  a second  order  effect, 

the  doping  is  thought  to  be  negligible  compared  to  the  doping 
31 

from  P (see  Ref.  1).  This  reaction  is  important,  however,  as 
the  primary  source  of  radioactivity  for  times  which  are  long 
compared  to  the  half-life  of  reaction  (3). 

From  the  data  given  in  Table  1,  we  can  calculate  the  prod- 
uction rate  of  the  desired  dopant.  The  total  number  of  Si  atoms/cm3 
is  found  to  be 

_ 2.33  (6.023  x 1022)  _ , „ n„22  „^„„y„„3 

N (Si ) — 25  0*6 — 4.99u  x 10  atoms/cm  . 

30 

The  density  of  Si  isotopes  is  then  found  to  be 

N (30Si)  = 0.0305  N(Si)  = 1.524  x 1021  atoms/cm3. 

The  theoretical  production  rate  for  phosphorus  is  then 
found  from  the  expression 


5 


d r 3 1 P 1 _ .,,30,,.,  _ r 30,.  • , . 31„., 

= N(  Si)  a [ Si  (n,  y)  Si] 

= 1.676  x 10  4 3^P  atoms/cm'Vn/cm^ 

O 

where  $ is  the  irradiation  fluence  of  thermal  neutrons  (n/cm  ) . 
Table  2 lists  the  theoretical  production  of  donors  in  Si  in  several 
useful  sets  of  units. 


TABLE  2. 


31 


P- Theoretical 


Production  Rates. 


Production  Rate  (K) 

Units 

-L 

1.676  x 10 

it 

1.676  x 10 

3.355 

3 2 

P atoms/cm  per  n/cm 

. 3 , -18  , 2 

P atoms/cm  per  10  n/cm 

^ ,nl8  , 2 

ppb  per  10  n/cm 



Knowing  the  concentration  of  donors  added  per  unit  fluence, 

it  is  now  possible  to  calculate  the  final  resistivity  after 

irradiation  and  annealing  from  the  starting  resistivity  if  the 

mobility  of  the  majority  carrier  is  known.  To  this  end  an  analysis 

of  mobility  data  for  Si  has  been  performed  on  data  taken  from 

2 3 

two  review  papers.  * Since  these  authors  quote  the  same  ref- 
erences but  give  different  values,  we  have  averaged  when  this  has 
occurred.  We  have  also  discarded  three  values  which  differed 

abnormally  from  the  mean,  p^n  = 1345,  ypn  = 1590,  and  1^=356 
2 

cm  /V-sec. 

The  result  of  the  average  of  all  other  data  is  summarized 
in  Table  3. 


r 


TABLE  3.  Typical  Experimental  Mobility  Values  in  Si 


= (U9^  +11)  cm  /V-see 
Dp  — 


y_  = (1396  + 6U)  cm  /V-sec 
Dn  — 


r = Ujj/vip  = 0.77  + 0.06 


r = 1.3  + 0.1 

n — 


u/u  = 2.826 

n p 


= Ceu  ]- 


16  -3 


[eynI 


(1.265  + 0.29)  x 10  cm 


= (253  + 5.7)  PPb 


(U.l+77  + 0.197)  x 1015  cm-3 


= (89.5^  ± 3.9h)  PPb 


These  expressions  will  be  useful  expressions  for  estimating  the 
ppb  added  by  neutron  irradiation  in  the  sections  to  follow. 

If  the  initial  and  final  carrier  concentrations  are  large  at 
room  temperature  compared  to  the  number  of  intrinsic  carriers, 
then  the  following  expressions  for  the  carrier  densities  in  terms 
of  donor  and  acceptor  concentrations  are 


p-type  starting  material 


P„  = na  ' Nn  . 


no  = ND  " NA’ 

O 


NA  ' NDf  ’ 


n£  " ND. 


n-type  starting  material 


p-type  ending  material 


n-type  ending  material 


then, 


n£  - n = Nn  - Np.  , 
£ O Dp  D * 
r o 


n-type  to  n-type 


nf  + p 

f 


NDf  ’ ND  » 
I o 


p-type  to  n-type 


V ' V 

f o 


p-type  to  p-type. 


But  Nn  - Nn  is  just  the  concentration  of  donors  added  by  trans- 
uf  uo 

mutation  which  we  will  call  Np  = Ki>  where  K is  the  production 

2 

rate  given  in  Table  2 and  $ is  the  fluence  (number  of  neutrons/cm  ) 
Therefore,  we  see  that  in  all  cases  for  which  n^  can  be  neglected, 
the  concentration  of  phosphorus  added  can  be  determined  by  the 
concentrations  of  uncompensated  carriers  and  the  material  type. 

We  can  determine  the  concentration,  Np,  of  phosphorus  atoms 
added  either  from  the  fluence  $ or  the  initial  and  final  resis- 
tivities and  types  using  the  following  expressions, 

Nd  = K$  = Hs.  355  x 10'18  71  $ (5) 

F n/ cni 


8 


Np  = l 


n 


n. 


n 


n 


K K 

-2.  + -E 


n. 


K 

_E 

PP„ 


K 

_E 

PPf 


n-type  -*■  n-type 


p-type  -*■  n-type 


p-type  -*■  p-type 


(6) 


where  = 89.54  ppb  ft-cm  and  = 253  ppb  ft-cm.  The  conversion 
between  ppb  and  atoms/cm3  is  1 ppb  = 4.996  x 1013  atoms/cm3. 

By  equating  the  two  expressions  above,  the  fluence  necessary 
to  obtain  a desired  final  resistivity  can  be  determined  from  the 
initial  resistivity  and  type  as 
K 


$ = < 


n t -1 
r~.  (pnf 


n. 


L) 


n p 


-1 


n. 


V 


-1 


I JV 


Cp, 


■l  -l, 

- ) 


Pf 


n-type  •+  n-type 

p-type  n-type 

p-type  -*■  p-type 


(7) 


The  resistivity  after  irradiation  and  annealing  as  a function 
of  fluence  can  also  be  calculated  for  those  cases  where  n^  is  not 
negligible  compared  to  the  majority  carrier  concentration. 

The  resistivity  vs  fluence  has  been  calculated  for  transmu- 
tation doped  samples  using  the  following  assumptions: 

1.  1004  electrical  activity  is  obtained  by  annealing. 

2.  The  mobilities  remain  constant  and  are  taken  as 

2 2 

Pn  = 1500  cm  /V-sec,  = 600  cm  /V-sec  in  the  following 
figures.  [S.M.  Sze , Physics  of  Semiconductor  Devices, 

Wiley,  New  York  (1969)]. 
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3.  The  shallow  donors  and  acceptors  are  completely  ionized 
at  room  temperature. 

4.  The  production  rate  of  phosphorus  is  given  by  K = 3.355 

18  2 

ppb/10  n/cm  , which  is  obtained  from  the  cross-section. 

5.  The  number  of  intrinsic  carriers  at  room  temperature  is 

in  _ x 

n^  = 1.391  x 10  cm  . [H.Y.  Fan,  Solid  State  Physics  1, 
283  (1955)]. 


From  charge  conservation,  the  complete  ionization  assumption, 
and  the  relation 

2 


np  = n. 


we  obtain 


P = i/2 [na  - nd  ♦ /vl/Tv], 

For  a production  rate  of  donors  given  by  K,  this  equation  becomes 
for  p-type  material 

(Na  - Nd  - K4>)  + - Nd  - KdO2  + 4n.*J  (8) 


p = 1/2 


If  we  define  tKe  degree  of  initial  compensation  kQ  by  the 


expression 


k = ND 
A 


then  the  degree  of  compensation  as  a function  of  fluence  is  given 
by 

k - ND  + K*  . 

For  complete  compensation,  k = 1;  for  no  compensation  k = o.  We 
define  one  more  quantity,  a,  the  degree  of  uncompensation  by 

a = 1-k  - NA  ' ND  ‘ K<t>  • p . ( 

- 5 ^ NA 
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p 


For  no  irradiation,  we  assume  that  pQ  >>n^  so  that 


a = 1-k 
o o 


NA  - ND  . Po 


(9b) 


N, 


N, 


A "A 

It  can  be  shown  by  substitution  of  Eq.  (9a)  into  Eq . (8)  that 


NA 

P = T- 


, m 2 , I 2n. 
a +.  / a + i 


2 -i 


N, 


(10) 


We  now  define  the  fluence  necessary  for  complete  compensation 
. 1 
Eq.  (9a) 


as  d>  . This  is  obtained  when  k = 1 or  a = o and  is  obtained  from 
c 


o = 1 - 1 = NA  - ND  - K<Jic 


TT 


or 


0 = N.-Nri  = p 

vc  A D *o 

K T~ 


It  can  now  be  shown  by  substituting  K 
last  equation  into  Eq.  (9a)  that 


(Na  * Nd)/^c  from  the 


a = 


(■■%)  (“(.) 


' a-  C1  • l). 


Then,  from  Eq.  (10) 
P 


(10)  r 

Vo  [(l  - | ) * / 


(1  - i ) ♦ 


(a  ■] 


But  from  Eq.  (9a)  with  <J> 

N . a = p 
A o *o 

therefore,  we  obtain 


= o 


11 


(11) 


It  can  similarly  be  shown  for  overcompensation  (n-type 
produced)  

[-(1-fe)+  v/tl ' !c)2  * U2) 

It  should  be  noted  that  the  number  of  carriers  is  symmetrical 
with  respect  to  the  fluence  on  either  side  of  if  * <f>  . It  should 
also  be  noted  that  the  number  of  free  carriers  as  a function  of 
fluence  is  independent  of  the  degree  of  initial  compensation. 

It  can  be  shown  that  Eqs.  (11)  and  (12)  form  a universal  doping 
curve  since  doubling  the  initial  concentration  doubles  the  nor- 
malized fluence  coordinate  (1  - ) , i.e.,  if  pQ  = 2 p ' , then 


P 


/k1  -fpw^y] 

C \ *0  / 


An  expression  for  the  resistivity  as  a function  of  fluence 
can  be  obtained  from  the  expressions 

1 


P = 


e(nyn  + pyp) 


and 


b = un/bp 


np  = n. 


12 


(13) 


By  substituting  n or  p as  a function  of  fluence  from  Eqs.  (11)  and 
(12)  the  resistivity,  p,  can  be  obtained  as  a function  of  fluence. 
This  expression  is  plotted  in  Figures  (1)  and  (2)  for  typical  p- 
type  float  zone  material  [the  parameters  are  given  on  Figures  (1) 
and  (2)  and  differ  from  those  on  Figures  (3)  - (6)]. 

The  sensitivity  of  resistivity  to  fluence  should  be  noted  near 
complete  compensation.  Furthermore,  in  the  expanded  scale  in  Figure 
2,  note  the  maximum  resistivity  occurs  at  a fluence  below  that  neces- 
sary for  100%  compensation.  This  is  a result  of  the  fact  that  the 
carrier  mobilities  are  not  equal.  The  intrinsic  resistivity,  , 
is  the  resistivity  obtained  at  exact  compensation. 

By  differentiation  of  Eq.  (13),  it  can  be  shown  that 

(14) 

where  pm  is  the  highest  resistivity  obtainable.  Also,  by  using 
a = neyn  + pepp  and  setting  do/dp  = o,  one  obtains  the  hole  concen- 
tration pm  = n^  )/ Vn/Pp  corresponding  to  the  maximum  resistivity  pm« 
It  can  be  similarly  shown  that  the  fluence,  <t>m,  to  obtain 
maximum  resistivity,  is  related  to  the  critical  fluence  to  obtain 
100%  compensation  by  using  Eq.  (11)  with  <{>  = 4>m  and  the  expressions 

p = p = n.  1 /p  /u  and  $ = p /K  to  give  the  expression 

v Km  i ¥ n p c 1 o 
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Resistivity  (s)  - cm) 


(1  - ♦/♦c) 


♦0.003  +0.006  +0.004  +0.002 


4.12  4.14  4.16  X 1016 


Fluence  (n/cm^) 

Figure  2.  Expanded  fluence  scale  plot  of  Figure  1 
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r 


= ' 1)(ni/K)/^n^ 


(15) 


which  depends  only  on  material  parameters  n^,  the  mobili-ty  ratio, 

and  the  phosphorus  production  rate  K.  This  difference,  therefore, 

* 

represents  a universal  constant  for  a particular  type  of  material. 
It  is  interesting  to  note  that  since  n^  and  K can  be  measured 
independently,  an  experiment  to  generate  Figure  2 gives  a new 
technique  for  measuring  the  drift  mobility  ratios  independent  of 
any  direct  mobility  measurement. 

This  model  can  be  extended  to  include  mixed  conduction  effects 
due  to  compensation  for  the  Hall  factor,  R^,  the  Hall  mobility, 

Upj,  and  the  apparent  number  of  free  carriers,  +_  1/R^e. 


■ 

From  equations 

(11)  and 

(12) , we  find  that 

■ Po 

- i ) + 

J f1  - 1 ) 2 * i 

f^M2] 

1 

L-J 

' T L 

*c 

V V *c  / 

l p0 1 J 

(16) 


where  (+)  is  for  p>n  and  (-)  is  for  n>p.  This  equation  represents 
the  actual  number  of  carriers  as  a function  of  fluence,  (f»  and 
the  critical  fluence,  $ , required  for  exact  compensation. 

Because  of  mixed  conduction  effects  near  exact  compensation, 

the  Hall  coefficient,  R^  is  suppressed.  A simple  extension  of 

4 2 

conventional  theory  after  using  the  relation  np  = n^^  gives 
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w 


r , 2 
P(P 


H 


2,  2, 

ni  b ) , 


rr:  2772 

e (p  + ni  b) 
n(ni2  - n2b2) 

L e(ni2  + n2b) 2 


p>n 


n>p 


(17) 


where  b = Pn/Up  is  the  ratio  of  drift  mobilities, 
relation  between  Hall  mobility  and  R^, 


Using  the  usual 


= — RH 


(+)  Rh>o 


and  equation  (13) , 


(-)  R <o 
H 


C 


epp(bni  + P2) 


n 

T~. 27 


ey  (bn  + ni  ) 


bn 

T“ ? 


for  R^>o 


Pne(bn  + n^  ) for  R^<o, 


(13) 


we  obtain 


( 


H 


f 2 2,  2. 

Up(P  - ni  b ) 

(p2  + n^b) 

Pn(n2b2  - n.2) 
L b (n2b  + ni2) 


for  R^>o 


and  for  R^<o. 


(18) 


Note  that  near  complete  compensation,  the  Hall  mobility  is 
considerably  suppressed  as  p + n^b  or  n + n^/b. 

The  apparent  free  carrier  concentration  obtained  from  the 
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Hall  coefficient  is,  near  compensation,  larger  than  the  actual 
carrier  concentration  since 


, 2 2,, 2 
(P  + ni  b) 


(19) 


2 2 2 
(n^b  + n.V 


These  equations  are  now  in  a form  from  which  y^,  R^,  and  n^ 
or  p^  can  be  obtained  as  a function  of  fluence  by  direct  sub- 
stitution of  Eq.  (16)  into  Eqs.  (17),  (18),  and  (19).  Programs 
have  been  written  to  generate  these  variables  as  a function  of 
fluence.  Data  calculated  from  the  code  is  shown  for  a typical 


case  [p  and  <(>  different  than  for  Figs.  (1)  and  (2)]  in  Figs. 

(3),  (4),  (5),  and  (6). 

Note  that  as  p + n^b;  -*■  o,  y^  -*•  o,  pH  -*■  °°.  This  is  not 

the  same  condition  as  obtained  for  the  resistivity  which  reaches 

,U 

a maximum  at  p = b n^. 

Although  the  above  equations  relate  the  measurable  room 
temperature  electrical  properties,  p,  y^,  n^,  and  p^  to  the 
fluence,  <f>,  no  direct  information  is  obtained  from  room  temperature 
electrical  measurements  on  the  ratio  of  boron  concentration  to 
phosphorus  concentration.  We  have  shown  that  the  free  carrier  con- 
centration can  be  compensated  in  principle  to  any  degree  of  pre- 
cision, however,  additional  information  from  low  temperature 
measurements  is  required  to  determine  if  the  desired  ratio  of 
boron  to  phosphorus  has  been  obtained.  The  question  remains, 
given  a particular  value  of  irradiation  fluence  ratio, 
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Fluence  (n/cirr) 


Figure  3 


Resistivity  vs.  fluence 


Fluence  (n/cnr) 


Figure  4.  Hall  mobility  vs.  fluence. 


Fluence  (n/cnr) 


Figure  6 


Actual  carrier  concentration,  p,  n,  an 
Hall  coefficient  carrier  concentration 
Pu*  nu*  vs»  fluence. 


how  closely  does  this  ratio  represent  the  boron  to  phosphorus 
ratio?  The  answer  is  easily  obtained.  The  conclusion  from  the 
analysis  to  be  presented  is  that  is  exactly  equal  to  the 

compensation  ratio  only  at  exact  compensation,  or  when  <J>  = <J>  . 

For  other  fluences,  this  ratio  is  not  exactly  equal  to  the 
compensation  ratio  but  is  a good  approximation  for  all  practical 
cases  of  interest.  To  demonstrate  this,  we  will  consider  two 
cases,  the  case  of  undercompensation  and  the  case  of  overcompensa- 
tion. 

Case  I:  Overcompensation,  <J>>4>  , p-type  to  n-type 

The  initial  compensation  ratio  in  the  p-type  material  before 
irradiation  is  defined  conventionally  as 

ko  ■ VNA 

while  after  irradiation,  the  compensation  ratio  is  defined  con- 
ventionally as 

kf  ' VNDf  • 

Since  for  all  starting  material  available  that  has  been  adequately 
zone  refined,  p >>n., 

p = Nfl  - Nn  = K<b 


where  again,  K is  the  production  rate  of  donors  by  the  transmutation 

process.  Also 

N - Nd  ♦ K cf> , 
t o 

therefore,  the  final  compensation  ratio  is  given  by 
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/ 


N, 


kf  = 


N 


D , 


\ + 

ND  + K<f> 


£ ,~0' 

By  rearrangement  and  using  the  relation  K4>c 
that 


N.-Nn  it  can  shown 
A o 


CVnd  - 


-1 


kf 


CVV  l)'1* 

($/$c) 

definition  of  k 

o’ 

(v1  - O’1 

+ 1 

(v1-  r 

♦ £ 

From  the  above  equation,  we  find  that 

v1  ■ ko  * ! u-m-  *>  v 

C 


(20) 


Note  that  when  <J>  = 4>  , kf  * 1.  Exact  compensation  can  be  deter- 
mined without  an  explicit  knowledge  of  the  starting  material 
compensation  ratio,  k , if  the  initial  carrier  conceptrat ion , pQ , 
and  the  phosphorus  production  rate,  K,  are  known.  In  general, 

v1  ■ vna  * 

Case  II:  Undercompensation,  cp  < 4>c , p-type  to  p-type 

In  this  case,  the  initial  and  final  compensation  ratios  are 
defined  conventionally  by  the  expressions 

kf  ■ vna-  k° ' vn*  • 

F o 


Therefore , 


Nd  + K* 

V s 0 

Kf  N 


°o  * K*c 


and  by  an  analysis  similar  to  Case  I,  we  obtain 
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(21) 


kf  = kQ  + (1  - kQ)  <j>/ <t>c ; <()CH  . (21) 

Again,  we  find  that  for  <J>  = <pc,  k^  » 1 but  for  other  ratios  of 

<P/<PC,  kf  f 4>/4>c - 

The  results  of  Eqs.  (20)  and  (21)  are  represented  graphically 

in  Figure  7.  In  this  figure,  the  final  ratio  of  phosphorus  to 

boron  concentrations,  /N^,  is  plotted  as  a function  of  the 

starting  material  compensation  ratio,  Nn  / N.,  for  different 

o 

fluence  ratios  4> /$  . It  is  seen  that  for  initial  compensation 


ratios  less  than  a few  percent 

* . \ 

IT 


(22) 


c A 

to  a very  good  approximation.  Starting  compensation  ratios  of  this 
order  of  magnitude  should  exist  for  multipass  zone  refined  Si 
because  of  the  differences  in  segregation  coefficients  for  boron 
and  phosphorus.  Even  when  the  initial  compensation  ratio  is  as 
high  as  20%,  the  error  in  estimating  the  final  compensation  ratio 
using  Eq.  (22)  is  no  more  than  20%  for  4>/<t>c  ^ 2,  an  overcompensa- 
tion of  up  to  100%.  Estimates  of  the  compensation  ratio  by  any 
known  experimental  technique,  including  low  temperature  electrical 
measurements  are  not  likely  to  exceed,  for  the  low  impurity  con- 
centrations with  which  we  are  here  dealing,  the  accuracy  which  can 
be  obtained  from  simple  room  temperature  electrical  measurements 
coupled  with  information  about  the  number  of  zone  passes. 

It  should  be  noted  in  concluding  this  section  that  the 
fluctuations  in  carrier  concentration  have  not  been  taken  into 
account  in  the  previous  equations.  This  will  be  discussed  further 


25 


UNDERCOMPENSATION:  ®/4>.  < I 

% 

k'un  = NT  = ko+  (l'ko)  *'*« 

OVERCOMPENSATION-  > I 


Figure  7.  Final  compensation  ratio  as  a function 
of  initial  compensation  ratio  for 
various  fluence  ratios,  <f>/<f>  . 


in  a later  section.  The  equations  presented  in  this  section  are 
exact;  however,  they  can  only  apply  to  local  concentrations  in 
the  wafer,  not  averages.  This  is  illustrated  by  considering,  for 
example,  our  definition  of  critical  fluence,  4>c.  Since  the 
initial  carrier  concentration,  pQ , is  a function  of  position  on 
the  wafer, 

P0  = po(x*  y)’ 

and,  from  our  previous  definition,  4>c  = pQ/K, 

<P c = <f>c  (x,  y). 

We  are  faced  with  the  perplexing  problem  that  different  parts 
of  the  wafer  require  different  fluences  to  reach  exact  compensation! 
If  we  irradiate  so  as  to  compensate  the  average  concentration 
exactly,  parts  of  the  wafer  will  be  undercompensated  and  other 
parts  will  be  overcompensated.  Under  such  circumstances,  the 
definition  of  exact  compensation  fluence,  $ , must  be  the  subject 
of  considerable  investigation.  We  have,  indeed,  identified  at 
least  a half  dozen  possible  definitions  for  the  critical  fluence. 

All  of  these  critical  fluences  differ  slightly;  however,  it  can 
be  shown  that  most  fall  within  +_  0.5%  of  each  other  for  typical 
starting  material  fluctuations. 

To  illustrate  some  of  the  problems  involved,  we  make  the 

hyphothesis  that  we  can  use  average  values  for  the  electrical 

measurements  to  determine  d>  . Then 

c 


for  p >>  n. 
*o  1 


or,  in  terms  of  average  resistivity, 
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c K p 

But  it  can  easily  be  shown  that  these  two  quantities  are  not 
equal  for  a working  definition  of  pQ  and  pQ  given  by 


Pq  max  + Po  min 


p + p 

Ko  max  Ko  mm 


Since  p = K /p  . then 
r o p o 


K.  -l 

p0  = “z2  (sr 


— + — ) 

0 p . ' 

Ko  max  o min 


= K p / K 

P 2 * P 

p p -ly 

o max  o mm  o 

It  can  be  shown  with  some  difficulty  that  the  error  intro- 
duced by  using  equation  (23)  to  define  <|>  compared  to  a different 
definition  for  the  critical  fluence,  $ ' , which  produces  equal 
numbers  of  carrier  fluctuations  on  either  side  of  intrinsic  can  be 


estimated  by  the  relation 
^c  1 I"  po  max  + 

^c  ^ L Po  min 


o mm 
5o  max 


Based  on  the  minimum  and  maximum  resistivity  values  for  samples 
VZ-083,  we  find  that 
<J>  - 

= 1.0025. 

*c 

A total  fluence  error  of  about  0.25%  is,  therefore,  introduced  by 


using  Eq.  (23)  and  the  above  definition  for  average  resistivity. 
The  definition  of  p given  above  is  also  suspect  since  the 


1 


true  average  (which  may  or  may  not  be  equal  to  estimates  made 
from  4-point  probe  data  in  a particular  experiment)  is 


p(x,  y)  = A'1  JJ  p(x,  y)  dxdy 
where  A = JJ"  dxdy.  The  previous  definition  for  p", 

p + p . 

— _ Kmax mm 


assumes  implicitly  a symmetrical  distribution  of  data  points  about 
the  mean.  Since  resistivity  is  not  directly  proportional  to 
carrier  concentration,  and  since  fluctuations  in  carrier  concen- 
tration depend  on  pull  rate,  rotation  rate  and  temperature 
fluctuations  during  growth,  it  is  very  unlikely  that  the  simple 
definition  for  p”  used  above  will  hold  exactly.  We  can,  however, 
be  reasonably  encouraged  from  the  estimate  of  critical  fluence 
error  introduced  by  these  approximations.  It  is  important,  how- 
ever, to  maintain  a certain  degree  of  caution  when  working  in 
the  regime  of  exact  compensation,  a regime  accessible  for  the 
first  time  only  by  transmutation  doping. 

III.  EXPERIMENTAL  TECHNIQUES  AND  APPARATUS 


In  this  section,  a description  of  the  apparatus  and  techniques 
required  to  irradiate  the  Si  and  to  evaluate  the  material  before 
and  after  irradiation  is  described. 

A.  Reactor  Description 

The  reactor  at  the  University  of  Missouri  is  a high  flux, 
light  water  moderated  research  reactor  designed  to  operate 
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on  highly  enriched  fuel.  It  has  a beryllium  reflector,  a 
graphite  moderator  and  forced  cooling  system  for  pool  and 
core.  The  reactor  is  operated  at  a power  level  of  10  MW. 

The  reactor  runs  routinely  at  10  MW  for  120  hours  per 
week,  52  weeks  per  year.  A comparison  of  the  peak  flux  of 
the  MURR  reactor  with  other  research  reactors  in  the  United 
States  is  shown  in  Figure  8. 

A cross-sectional  view  of  the  reactor  core  and  various 
experimental  stations  is  shown  in  Figure  9.  One  of  the 
major  advantages  of  the  MURR  reactor  for  Si  irradiation  is 
the  easy  access  to  the  sample  positions  in  the  pool  and 
reflector. 

At  present,  two  2"  dia.  x 30"  long  silicon  sample 

irradiation  positions  with  sample  rotators  are  located  in  the 

graphite  reflector  in  positions  1-2  and  H-2  in  Figure  9.  The 

13  2 

average  flux  at  these  positions  is  about  3 x 10  n/cm  /sec. 

The  pneumatic  tubes  used  for  activation  analysis  rapid 

sampling  handling  are  also  shown  in  Figure  9. 

Beyond  the  graphite  positions  is  the  S-basket  with  seven 

3"  dia.  Si  holes  which  can  be  rotated  manually  or  mechanically. 

13  2 

The  flux  in  the  S-basket  varies  from  10  n/cm  /sec  to 

ml2  / 2. 

10  n/cm  /sec. 

For  the  very  low  fluences  required  for  the  exact  com- 
pensation of  detector  Si,  a variable  flux  pool  facility  was 
developed  and  constructed.  The  samples  in  this  position  are 
also  mechanically  rotated.  The  flux  in  this  facility  is  a 
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Figure  8.  Comparison  of  research  reactor  peak 
flux  values  in  the  United  States. 


Figure  9. 


Flux  trap,  beamports,  reflector  and 
bulk  pool  facilities. 
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maximum  in  the  present  close-in  position  and  is  about 
11  2 

7 x 10  n/cm  /sec.  This  flux  can  be  reduced  considerably 
by  moving  the  facility  further  out  in  the  pool.  The  approx- 
imate position  of  the  bulk  pool  Si  facility  is  shown  in 
Figure  9. 

A typical  flux  profile  for  the  1-2  position  is  shown 
in  Figure  10.  Silicon  ingots  up  to  10"  length  are  positioned 
in  the  linear  portion  of  this  flux  and  flipped  end  for  end 
at  precisely  50%  of  the  total  fluence.  The  flux  profiles  in 
other  positions  are  approximately  the  same  shape  but 
tend  to  flatten  at  the  peak  at  positions  further  from  the 
core.  In  the  variable  flux  facility,  the  peak  is  suffi- 
ciently flat  to  make  the  peak  of  the  flux  profile  useful 
over  a 5"  length. 

B.  Bulk  Pool  Variable  Flux  Irradiation  Facility 

A calculation  of  the  fluence  required  to  compensate  a 

1 2 

typical  residual  boron  concentration  of  about  3 x 10  atoms/ 

3 16  2 

cm  yields  a fluence  of  1 to  2 x 10  n/cm  . A typical  loading 

or  unloading  time  of  five  minutes  and  a fluence  accuracy  of 

0.5%  restricts  the  irradiation  flux  to  a value  between 
11  12  2 

1 x 10  to  1 x 10  n/cm  /sec.  Since  this  flux  is  an  order 
of  magnitude  lower  then  that  obtained  in  the  S-Basket.  a 
special  irradiation  facility  w«s  constructed  for  this  project. 
The  general  location  of  this  "Bulk  Pool  Variable  Flux  Facility" 
is  shown  in  Figure  9. 

Figure  11  shows  the  construction  details  of  the  movable  box 
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which,  is  held  onto  the  S-Basket  by  two  rigid  rods  which 
pass  through  holes  A.  (The  sample  holding  tube  is  oriented 
in  the  pool  so  that  "Top"  in  Fig.  11  points  toward  the  top 
of  the  pool).  Fig.  12  shows  a detail  of  the  support  rods 
while  Fig.  13  shows  the  detector  holder  which  inserts  into 
the  main  cylinder  in  Fig.  11.  The  sample  holder  is  a cylinder 
which  rotates  inside  the  detector  hole  in  Fig.  13.  It  should 
be  noted  that  the  neutron  detector  spirals  concentrically 
around  the  sample  position  so  that  a very  accurate  representa- 
tion of  the  flux  on  the  sample  is  continuously  monitored.  The 
thermal  to  fast  flux  in  this  facility  is  estimated  to  be 
better  than  30:1. 

C.  Neutron  Flux  Detection  and  Integration  System 

The  detector  discussed  in  the  last  section  is  a self 
powered  detector  of  the  thin  wire  design. 

103 

Thermal  neutrons  are  captured  in  the  thin  Rh  wire 
(1001  abundance) . A nuclear  reaction  is  induced  which  is 
proportional  to  the  flux  of  the  type 


103D.  , s 104  , 104  , x - - 

Rh(n,  y)  Rh  -*■  Pd  + 6 . 

The  half  life  for  the  3 emission  is  4.2  sec.,  therefore, 
this  detector  is  able  to  follow  the  rather  slow  flux 
variations  caused  by  movement  of  the  reactor  control  rods. 

The  electrons  from  this  reaction  are  collected  on  a con- 
centric shield  to  produce  a current  proportional  to  the  flux. 
This  current  is  then  integrated  and  an  alarm  set  off  at  a 
predetermined  preset. 
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Lgure  12.  Support  rods  for  variable  flux  facility 


Figure  13.  Detector  holder  for  variable  flux  facility. 
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Considerable  effort  has  been  spent  calibrating  this 
system.  Flux  wires  of  Co  doped  A1  have  been  used  to  measure 
the  flux  by  the  radioactivity  generated  by  a given  fluence 
and  compared  to  the  observed  counts  from  the  integrator  system. 
Two  types  of  flux  wire  have  been  used,  N.B.S.  and  Cohn. 

Although  considerable  variation  due  to  flux  wire  counting 
statistics  and  Co  concentration  are  observed,  the  running 
mean  of  the  two  sets  of  flux  wires  agree  to  an  accuracy 
approaching  that  of  the  integration  system  accuracy  (+T).2%). 
This  is  illustrated  in  Table  4. 


I 


TABLE  4.  Integrator  System  Calibration  Data. 

Irrad.  # Wire  Fluence /Count  /^Deviation  from  Running 

(n/cm2 /count)  Last  Mean  Mean 


Cohn 

4 

33-34 

1.557  x 1012 

- 0.128? 

12 

1.557  x 10 

8 

33-34 

1.610 

+ 3.27 

1.594 

9 

33-34 

1.568 

+ 0.577 

1.578 

11 

33- 34B 

1.549 

- 0.641 

1.576 

12 

33-34 

1.511 

- 3.079 

1.559 

5 

NBS 

1.553  x 1012 

- 0.830? 

12 

1.553  x 10 

8 

NBS 

1.482 

- 5.364 

1.518 

11 

NBS 

1.664 

+ 6.258 

1.5  66 

100  x NBS  - Cohn 

(NBS  + Cohn ) / 2 

0.449? 

We  believe 

we  have  demonstrated 

that  the 

fluence  as 

determined  by  the  integration  system  is  accurately  known  and 


reproducible  to  about  +T).5%. 
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D. 


Neutron  Activation  Analysis  Systems 


Existing  neutron  activation  analysis  (NAA)  facilities 

were  employed  for  measuring  elemental  impurities  in  the 

silicon  ingot  starting  materials.  Two  NAA  irradiation 

facilities  are  available:  a pneumatic  tube  transfer  system 

into  a neutron  flux  of  about  8 x 10  n/cm  /sec  which  permits 

the  measurement  of  relatively  short-lived  activation  products 

and  an  outer  reflector  position  at  approximately  1 x 10^n/ 

2 

cm  /sec  for  long  term  irradiations.  The  gamma-ray  spectro- 
scopy instrumentation  consists  of  a lead  shielded  50  cc,high 
resolution  Ge(Li)  detector  coupled  to  a 4096  channel  MCA 
with  computer  compatible  magnetic  tape  output.  An  I.B.M. 

370  computer  is  used  for  data  reduction.  These  facilities 
provide  "state  of  the  art"  capability  and  reliability  for 
neutron  activation  analysis  of  trace  element  impurities  in 
silicon. 

E . Four-point  and  Thermal -Probes 

As  discussed  in  Section  II,  considerable  information 
can  be  obtained  about  the  nature  of  the  compensation  process 
by  measuring  resistivity  and  majority  carrier  type  after 
typical  majority  carrier  mobility  is  established  for  the 
ingot  or  wafer. 

Thermal  probe  measurements  of  majority  carrier  type  are 
made  in  a conventional  manner  using  a variable  temperature 
soldering  pencil  tip  for  the  hot  probe  and  a steel  tip  cold 
probe.  The  temperature  gradient  generated  signal  is  fed 
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into  a Keithley  Model  602  electrometer 5 operated  in  the 

center  zero  current  measuring  mode  usually  on  a range  of 
- 8 - 7 

10  to  10  A.  The  hot  probe  is  connected  to  the  electrometer 
low  while  the  cold  probe  is  connected  to  the  high  so  that  a 
positive  deflection  denotes  p-type  and  a negative  deflection 
denotes  n-type.  The  type  of  the  unknown  sample  is  also 
compared  with  two  standard  samples  obtained  from  a multipass 
zone  refined  ingot  (p-type)  and  a neutron  transmuted  and 
annealed  sample  (n-type).  Agreement  of  the  thermal -probe 
with  Hall  effect  measurements  has  been  excellent  up  to 
resistivities  approaching  100,000  Jl-cm  once  experience  is 
gained  with  the  standard  samples.  Beyond  this  resistivity, 
mixed  typing  is  often  encountered  suggesting  that  nonuniform 
doping  is  observed  before  the  disappearance  of  the  thermoprobe 
signal  is  observed. 

Several  four-point  probe  heads  of  different  probe  spacing 
and  probe  pressure  have  been  tried  with  excellent  results  in 
resistivity  ranges  usually  inaccessible  to  four-point  probe 
measurement.  The  system  consists  of  a Keithley  Model  225 
constant  current  source5,  a Hewlett  Packard  Model  3465A 
digital  multimeter^,  and  an  Alessi  Industries  ATP-1  Test 

7 

Probe  Fixture'  which  has  been  modified  to  nave  a 12"  vert- 
ical travel  and  a sample  translation  stage.  Two  types  of 
four-point  probe  heads  have  been  tried,  an  Alessi  Model 
A4P-25  with  a probe  spacing  of  0.025M,  a probe  load  of  70- 

O 

180  gm;  and  a Signatone  probe  with  a spacing  of  0.065",  a 
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tip  radius  of  0.1  mil  and  a probe  load  of  180  gms.  Both 

probe  heads  have  tungsten  carbide  tips.  Although  both 

probe  heads  work  well,  the  larger  probe  spacing  is  more 

forgiving  in  current  selection  and  care  of  measurement  and 

has  been  used  except  for  very  small  samples  where  edge 

effects  become  a problem.  The  usual  probe  correction  factors 

have  been  applied  to  obtain  the  resistivity  from  the  measured 
9-11 

resistance. 

To  account  for  nonuniform  probe  spacing,  the  probes 
were  lowered  onto  an  aluminized  microscope  slide  and  the 
spacing  measured  with  a traveling  microscope.  This  procedure 
was  repeated  several  times  and  a correction  term  generated 

q 

from  theory.  This  correction  factor  was  also  checked  by 
measuring  the  resistivity  of  a transmutation  doped  and  an- 
nealed sample  which  had  been  calibrated  by  the  Van  der  Pauw 
1 2 

technique.  Agreement  was  found  to  be  better  than  1%  for 
a 5 n-cm  sample.  All  samples  are  measured  in  a dark  box. 

Since  the  A.S.T.M.  Standard  for  resistivities  over 
2000  ft-cm  suggested  that  measurements  could  not  be  made  to 
an  accuracy  of  better  than  +_  15%,  extreme  care  was  taken 
to  check  the  probe  calibration  factor  periodically  for  high 

TP'i  i ^ t i VI  wirtocn-romont’C  ^ ^ 

iCJijtivicr  in  m a w>.>v<4  • 

An  ingot  which  had  been  cut  into  five  bulk  samples 
(p  s,  2500  fi-cm)  and  characterized  by  four-point  probes  at 
Rockwell  International  was  measured  at  MURR.  Care  was  taken 
to  use  a current  sufficiently  low  that  sample  heating  did  not 
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affect  the  measurements.  The  results  of  I -V  plots  for  two 
sides  of  sample  VZ-072-#S  are  shown  in  Figures  14  and  15.  The 
nonlinearity  in  these  measurements  at  high  currents  is  due 
to  sample  heating. 

Note  that  the  linear  portions  of  the  curves  do  not  go 
through  the  origin.  This  offset  voltage  can  be  measured 
on  the  voltmeter  in  the  absence  of  current  and  is  observed 
to  vary  with  probe  pressure,  surface  condition  of  the  sample 
and  position  of  the  probes  on  the  sample.  Although  this 
effect  is  negligible  for  5 Q-cm  material,  the  offset  voltage 
can  become  an  appreciable  part  of  the  signal  voltage  for  high 
resistivity  samples.  It  was  also  noticed  that  the  slope 
near  the  origin  is  independent  of  the  offset.  This  fact 
allows  us  to  make  correct  measurements  up  to  very  high 
resistivities  if  proper  data  handling  is  used.  The  problem 
is  overcome  by  measuring  both  the  positive  and  negative 
voltages  resulting  from  positive  and  negative  supply  currents 
to  determine  the  effective  resistance.  This  resistance  is 
then  equal  to  [V(+)  - V(-)]/[I(+)  - I(-)].  This  resistance 
determines  the  slope  of  the  lines  in  Figures  14  and  15  and 
the  product  of  this  slope  times  the  probe  factor  determines 
the  resistivity. 

Results  of  a series  of  measurements  on  the  five  bulk 
samples,  VZ-072-#l  through  #5,  are  compared  to  results  obtained 
by  Dr.  Arst  at  Rockwell  and  are  shown  in  Table  5. 

It  is  seen  that  the  agreement  between  the  two  laboratories 
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TABLE  5*  A comparsion  of  four  point  probe  resistivity  measurements  made 
at  Rockwell  and  MURR. 


Sample 

Rockwell  Data 
Resistivity  (ft-cm) 

MURR  Data 
Resistivity  (ft-cm) 

100  Pr-Pm 

pm 

0l(T) 

2.29  x 103 

2.27  x 103 

+0.88 

01(B) 

2.42 

2.1+3 

-0.1+1 

02(T) 

2.43 

2.1+5 

-0.82 

02(B) 

2.53 

2.53 

0.00 

03(T) 

2.50 

2.52 

-0.99 

03(B) 

2.57 

2.58 

-0.39 

0l+(T) 

2.57 

2.58 

-0.39 

04(b) 

2.55 

2.56 

-0.39 

05(T) 

2.55 

2.55 

0.00 

05(B) 

2.53 

2.1+9 

+1.61 

is  quite  good  for  this  resistivity  range  the  average 
deviation  being  about  +_  0.61.  The  agreement  between  the 
four-point  probe  measurements  has  been  observed  to  decrease 
with  further  increases  in  resistivity.  On  comparing  both 
laboratories'  four-point  probe  measurements  with  Van  der 
Pauw  measurements  on  samples  with  resistivities  in  the 
7000-8000  fi-cm  range  we  find  the  deviation  at  MURR  from  Van 
der  Pauw  to  be  +_  1.5%  while  the  deviation  of  the  Rockwell 
measurements  is  + 5%. 


This  increased  error  with  resistivity  seems  to  be 
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following  a power  law  dependence  given  by  the  expression 

100  = Kpn  . 

P 

Although  the  constant  K is  smaller  for  our  data  than 
for  that  given  by  the  A.S.T.M.,  the  value  for  n is  in 
reasonable  agreement.  We  find  the  following  to  hold  for  the 
two  sets  of  data: 


100 

Ap 

P 

= 0.1  P0-71 

(A.S.T.M.) 

100 

Ap 

P 

= 0.001  p0-81 

(MURR) . 

The  physical  meaning  of  these  equations  is  obscure, 
however,  a similar  power  law  is  observed  for  the  calibration 
curves  for  spreading  resistance  measurements}'"  The  difference 
between  the  A.S.T.M.  and  MURR  values  for  K is  not  understood, 
however,  the  A.S.T.M.  data  is  taken  for  both  n-  and  p-type 
Si  in  various  orientations  while  the  MURR  data  is  taken  on 
p-type  only  on  a [ill]  face.  The  upper  expression  is  useful 
for  estimating  the  maximum  expected  error  for  intrinsic 
material  at  230,000  ft-cm.  We  find  an  expected  error  of  about 
+_  221  at  this  resistivity  value.  This  estimated  error  is  in 
reasonable  agreement  with  deviations  in  resistivity  for 
samples  neutron  irradiated  to  the  same  fluence  but  before 
annealing  (nearly  intrinsic  but  slightly  p-type). 

F.  Minority  Carrier  Lifetime  Apparatus 

The  minority  carrier  lifetime  of  all  ingots  before 
irradiation  and  some  after  irradiation  and  annealing  has 
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been  measured  by  the  photoconductive  decay  technique.  Ingots 
are  coated  with  Aquadag  solution  on  each  end  and  held  in  a 
small  lab  vise  to  which  braided  copper  shielding  wire  has 
been  flattened  and  attached  to  the  vise  jaws  for  contacting. 
The  contacting  area  is  shielded  from  the  light  to  prevent 
photovoltaic  effects  from  the  contacts.  This  effect  can  be 
detected  if  observed  as  a linear  decay  of  the  junction  ca- 
pacitance of  a blocking  contact. 

A small  power  supply  provides  the  bias  voltage  across 
the  sample  and  a variable  series  load  resistor.  This  resistor 

is  adjusted  for  maximum  signal  before  data  is  taken,  an 

14 

essential  procedure  for  interpretation  of  the  decay. 

Optical  excitation  is  provided  by  a General  Radio  Strobotac 
model  1531-AB^  whose  pulse  time  is  short  compared  to  the 
photoconductive  decay  time  constant.  The  signal  is  observed 
with  a Tektronix  Type  547  oscilloscope  with  a dual  trace 
plug-m  unit. 

Typical  decay  voltages  as  a function  of  time  are  shown 

in  Fig.  16.  It  can  be  shown  that  the  voltage  signal  as  a 

function  of  time  for  the  condition  that  R(load)  = R(sample) 

. 17 

is  given  by 

(£ ) . *~t/T 

where  (An/nQ)o  is  the  maximum  change  in  minority  carrier 


VCt) 

bias 
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concentration  and  t is  the  carrier  lifetime.  For  times 


long  compared  to  the  optical  excitation  decay 

e~t/T 

-Vs1  -(s2-) 

bias  ' no  ' o 


Therefore,  a semilog  plot  of  V(t)/Vbias  has  a slope  of  t 
where  x is  the  minority  carrier  lifetime. 


Care  has  been  taken  not  to  optically  overmodulate  the 


samples  since  this  is  known  to  produce  exaggerated  lifetimes. 


The  procedure  we  use  is  to  reduce  the  light  intensity  until 


the  observed  decay  time  constant  becomes  a minimum.  It  is 


also  important  that  the  carrier  sweep  out  time  given  by 

_ T 2 


t = L‘ 


y^bias 


where  L is  the  sample  length  and  y the  carrier  mobility,  be 


long  compared  to  the  decay  time, 


Trapping  effects  have  been  detected  in  some  samples  by 


observing  a decrease  in  the  decay  time  when  the  sample  is 


siiriul  tanennsl  v illuminated  from  a dc  tungsten  light  source 


(ftw<E  ) . In  these  cases,  the  lifetime  is  assumed  to  be  that 


observed  during  dc  illumination  if  the  reduction  in  decay 


time  saturates. 


Circuitry  has  been  developed  to  eliminate  the  need  to 


take  data  from  the  exponential  photoconduct ive  decay  displayed 


on  the  oscilloscope.  A variable  time  constant,  variable 


pulse  height  exponential  generator  has  been  constructed  which 
provides  a decay  time  constant  whose  value  can  be  obtained 
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from  the  setting  of  a ten-turn  pot.  This  exponential 
pulse  is  triggered  from  the  strobotac,  inverted,  and  added 
to  the  photoconductive  decay  signal.  The  oscilloscope  then 
functions  as  a null  detector.  A straight  line  signal  in- 
dicates that  the  exponential  generator  has  matched  the  photo- 
conductive  signal  in  both  amplitude  and  time  constant. 
Deviations  from  non-exponential  decays  and  trapping  effects 
are  easily  detected  with  this  system.  A schematic  of  the 
system  is  shown  in  Figure  17  while  a circuit  diagram  for  the 
exponential  pulse  generator  is  shown  in  Figure  18. 

G.  High  Impedance  Van  der  Pauw  Resistivity  and  Hall  System 


(Electronics,  Vacuum  System,  and  Helitran) . 

Resistivity  and  mobility  can  be  determined  as  a function 

of  temperature  from  about  8°K  to  room  temperature  by  this 

1 8 

system.  Temperature  control  is  obtained  from  a Helitran 
and  temperature  controller.  The  system  is  capable  of 
measuring  resistances  of  high  impedance  samples  up  to  a 
resistivity  of  about  10'L4  fi-cm.  A dc  guarded  lead  approach-1^ 
in  which  each  sample  lead  guard  is  driven  by  a separate  high 
input  impedance  unity  gain  amplifier  (Keithley^  Model  602) 
to  the  voltage  on  the  inner  lead  reduces  leakage  currents 
and  improves  the  system's  time  constant.  The  circuitry  for 
this  approach  is  shown  in  Reference  19.  The  digital  volt- 
meter and  constant  current  source  from  the  four-point  probe 
system  are  also  used  in  this  system.  A switching  circuit 
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CARRIER  LIFETIME  MEASUREMENT  APPARATUS 


EXPONENTIAL  DECAY  GENERATOR 


electrometer  and  current  source  allows  either  conventional 

12 

resistivity  and  Hall  effect  or  Van  der  Pauw  measurements 
to  be  made.  This  circuit  is  shown  in  Figure  19. 

A vacuum  system  has  been  constructed  to  provide 
evacuation  for  the  Helitran  system  and  for  the  vacuum  an- 
nealing. The  system  uses  an  oil  diffusion  pump,  mechanical 

forepump  and  liquid  nitrogen  cold  trap.  This  system  provides 

_ 7 

a vacuum  of  5 x 10  torr. 

The  magnetic  field  for  Hall  effect  measurements  is 

20 

provided  by  a Varian  Model  V4004  four-inch  magnet  system. 

The  magnetic  field  has  been  calibrated  to  +_  1.3%  using  both 
a rotating  coil  Gauss  meter  which  has  been  compared  to  an 
NMR  probe  and  a calibrated  semiconductor  Hall  probe.  The 
semiconductor  Hall  probe  is  used  to  set  the  magnetic  field 
using  circuitry  which  provides  digital  readout  of  the  field. 
This  circuitry  is  shown  in  Figures  20  and  21.  The  magnetic 
field  for  Hall  effect  measurements  in  this  report  is  6.195  kG. 
H.  Photoconductivity  System 

Because  the  transmutation  doping  studies  have  been  per- 
formed primarily  on  undoped  material,  the  necessity  of  per- 
forming photoconductivity  studies  to  obtain  the  degree  of 
compensation  or  the  photoconductive  spectral  responses  have 
not  been  needed.  This  system  has  been  developed,  however, 
in  anticipation  of  further  IR  detector  spectral  response 
work.  The  system  consists  of  a glow  bar  light  source,  a 
Perkin-Elmer  double  pass  prism  monochrometer  and  electronics, 
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POWER 


DIGITAL  KGAUSS  METER 


two  thermocouple  detectors  and  exit  optics.  The  prism 
and  windows  are  KRS-S  material  (1-50  ym)  . The  exit  beam 
from  the  monochrometer  which  is  diverging  is  intercepted  by 
a spherical  mirror  to  produce  a parallel  beam  of  light.  Part 
of  this  parallel  beam  is  interrupted  by  a diagonal  plane 
mirror  and  focused  on  a reference  thermocouple  detector 
with  a KRS-5  window.  The  other  part  of  the  beam  is  focused 
through  a second  window  by  another  spherical  mirror  onto  the 
photodetector  sample  mounted  in  the  Helitran.  By  taking  ratios 
of  these  two  signals  the  spectral  response  of  the  photodetector 
can  be  obtained  normalized  to  the  light  source  and  system 
response. 

I . Miscellaneous:  Sample  Preparation  Equipment,  Annealing 


System,  Irradiation  of  Silicon  Post-Irradiation 
Handling 

Samples  are  cut  into  wafers  either  by  a diamond  wire  saw 
or  diamond  disk  cut-off  saw.  The  latter  has  been  found  more 
suitable  for  silicon.  Samples  are  heavily  etched  after 


cutting  to  remove  s u 1 1 a v~e  u a ma g c in  a conventional  HT/HNO^ 


i u n /mtr\  / 

* r/ 


acetic  etch.  Samples  are  then  ground  using  various  grades 
of  SiC  and  Al^O^  grits  and  finally  washed  in  deionized  water, 
acetone,  methanol,  and  trichloroethylene. 

Samples  are  annealed  in  a tube  furnace  inside  a Spectrasil 


quartz  tube.  The  sample  may  be  annealed  in  either  a vacuum  of 
•7 


5 x 10  torr.  or  in  a flowing  inert  gas  atmosphere.  Sample 
surfaces  are  protected  on  each  side  with  a wafer  of  float 
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zone  silicon  to  prevent  contamination  by  the  quartz.  The 
furnace  is  capable  of  temperatures  of  up  to  1000°C.  The 
temperature  is  regulated  to  +_  0.5 °C. 


Samples  to  be  irradiated  are  placed  between  silicon 
spacers;  the  whole  stack  of  samples  and  spacers  is  then 
wrapped  in  aluminum  foil.  The  flux  wires  are  now  wrapped 
around  the  foil  and  the  whole  package  is  placed  in  an 
aluminum  can  which  is  welded  shut  to  protect  the  silicon 
from  the  pool  water.  Radioactivity  is  monitored  after  a 
cool-off  period  in  the  reactor  pool.  For  the  very  light 
irradiations  in  this  project,  radioactivity  decays  to 
background  very  quickly  and  has  usually  been  attributed  to 
surface  contamination  where  encountered. 

A considerable  effort  has  been  spent  studying  surface 
contamination  after  irradiation  but  before  annealing.  By 
careful  handling,  surface  contamination  has  usually  not 
been  detected  before  annealing  except  for  a Au  contamination 
on  one  sample  from  sawing  in  another  laboratory  before  our 
saws  were  operational.  A few  samples  have  also  shown  traces 
of  Na  from  handling.  This  has  been  eliminated  by  a light 
HF  cleaning  after  irradiation  but  before  annealing.  We  are 
confident  that  our  samples  are  usually  extremely  clean 
prior  to  annealing  based  on  activation  analysis  and  electrical 
measurements  after  annealing. 

IV.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 

thlS  SeCUOn'  " Pri0r't0-irradiati°n  "aterial  eValUaliW"  I 
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irradiation  and  annealing  experiments  are  discussed.  We  have 
noted  several  differences  in  our  results  from  previously  published 
work  and  will  discuss  these  differences  in  the  following  sub- 
sections. A preliminary  activation  energy  for  the  thermal  ion- 
ization of  the  irradiation  produced  donors  has  been  measured. 

We  have  been  able  to  determine  that  the  best  single  annealing 
temperature  for  NTD-Si  is  at  above  800°C  for  one  hour  followed 
by  a very  slow  cool  down.  We  have  also  made  a preliminary 
assessment  of  the  effects  of  fast  to  thermal  neutron  ratios. 

Although  material  from  a variety  of  Si  producers  has  been 
evaluated,  we  have  found  that  control  by  compensation  on  the  p- 
type  side  of  type  conversion  from  5000  ft-cm  to  intrinsic  is 
impossible  in  the  material  we  have  investigated  thus  far.  We  will 
present  data  which  suggests  that  residual  oxygen  is  the  source  of 
difficulty. 

W'e  have  also  demonstrated  very  good  control  of  resistivity  on 
the  n-type  side  of  type  conversion.  We  have  shown  that  the  final 
resistivity  can  be  controlled  to  at  least  +_  2%  using  only  a 
preliminary  check  of  ingot  mobility  and  four-point  probe  measure- 
ments. We  have  sound  experimental  evidence  at  this  point  which 
suggests  that  it  may  be  possible  to  achieve  this  degree  of  control 
up  to  n-type  resistivities  as  high  as  100,000  n-cm.  This 

- 4 

corresponds  to  an  excess  donor  concentration  of  'v  8 x 10  ppb 
10  3 

or  about  4 x 10  donors/cm  . It  may  thus  be  possible  to 
compensate  residual  boron  concentrations  of  about  this  order  of 
magnitude  or  larger. 
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One  of  the  more  exciting  possibilities  to  develop  out  of 
the  present  work  is  residual  oxygen  decoration  by  damage  defects. 
Circumstantial  evidence  suggests  that  oxygen  detection  in  Si  can 
be  pushed  at  least  three  orders  of  magnitude  below  present  detection 
limits.  If  the  preliminary  results  can  be  substantiated  by  EPR 
and  optical  measurements,  considerable  progress  in  Si  purification 
can  be  made  as  a result  of  this  technique. 

A sound  experimental  approach  has  been  developed  for  deter- 
mining the  amount  of  boron  in  the  samples  prior  to  irradiation, 

which  is  available  to  those  IR  detector  manufacturers  who  perform 
their  own  Ga  or  In  doping.  We  feel  that  this  can  be  done  success- 
fully with  high  accuracy  only  before  the  optically  active  impurity 
is  added  to  the  ingot.  Although  other  techniques  are  available 
to  determine  boron  concentration  after  the  addition  of  Ga  or  In, 
none  rival  in  accuracy  the  technique  of  multizone  passes  and  four- 
point  probing  prior  to  optical  activation  by  Ga  or  In  doping. 

We  have  also  been  able  to  establish,  both  theoretically  and 
experimentally,  the  resistivity  uniformity  to  be  expected  as  a 
function  of  distance  across  the  wafer  for  various  compensation 
ratios.  This  is,  perhaps,  the  most  important  single  result  of 
the  present  work  since  results  clearly  show  that  compensation  by 
nuclear  transmutation  offers  a significant  and  considerable 
improvement  over  conventionally  compensated  material  in  the 
compensation  range  useful  for  IR  detectors. 

We  have  also  been  able  to  relate  the  final  resistivity 
f luctuatuions  observed  for  various  compensation  ratios  to  the 
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initial  boron  fluctuations  in  starting  material  in  a quantative 
way.  It  is,  therefore,  possible  for  the  first  time  to  place 
criteria  on  the  boron  fluctuations  as  determined  by  the  desired 
final  wafer  uniformity  and  compensation  ratio.  The  theoretical  and 
experimental  results  we  have  obtained  indicate  that  nuclear  trans- 
mutation compensation  will  be  a necessary  technique  for  uniformity 
criteria  now  envisioned  for  Si  IR  detector  array  integration. 

A.  Materials  Evaluation  Before  Irradiation 

All  the  material  investigated  has  been  multizone  re- 
fined float  zone  silicon  with  resistivities  from  2500  ft-cm 
up  to  25000  ft-cm.  The  material  is  p-type  and  not  inten- 
tionally doped.  The  dominant  acceptor  is  undoubtedly  boron. 

All  of  the  material  is  nominally  dislocation  free  and  pulled 
in  a [111]  direction  except  the  Rockwell  material  which  is 
dislocated  and  pulled  in  a [100]  direction.  Some  preliminary 
electrical  data  on  a Wacker's  [100]  pulled  Ga  doped  sample 
have  been  taken.  Neutron  reflection  rocking  curves  show 
reflection  line  widths  of  about  1 sec  of  arc  for  the  dis- 
location free  material  and  about  1 min  of  arc  for  the  dis- 
located material.  Material  investigated  has  come  from  Dr. 

Arst,  Rockwell;  Topsil,  Wackers,  Monsanto,  Hughes,  and  a Texas 
Instruments  Lopex  sample  supplied  by  Dr.  Spry  of  AFML.  We  have 
added  an  in-house  number  to  each  ingot  or  wafer  (MURR  #)  and 
have  cross-referenced  these  numbers  with  whatever  sample  des- 
ignations have  been  assigned  to  the  samples  by  the  suppliers. 

Much  effort  has  been  expended  in  trying  to  purchase 
float  zone  silicon  which  is  doped  with  either  In  or  Ga. 
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Approximately  thirty  manufacturers  of  silicon  were  contacted 
by  phone  and  requests  for  bids  sent  out  to  ten,  however,  no 
bids  were  received  for  this  type  of  material.  We  have,  how- 
ever, received  free,  doped  samples  from  a number  of  sources 
and  greatfully  acknowledge  their  help.  We  have  received 
Czochralski  grown  Si:Ga  doped  to  a variety  of  levels  in  the 
range  of  10^  to  10^  cm  ^ from  Drs.  M.  Bressler  and  C.  M. 

Parry,  Aerojet  ElectroSystems  Co.  We  have  also  received  a 
supply  of  Wacker's  Ga  doped  vacuum  float  zone  material 
from  Dr.  M.  Schlacter,  Aeronutronic  Ford.  Dr.  Arst  of 
Rockwell  has  also  supplied  us  with  Ga  doped  float  zone. 

Table  6 summarizes  sample  characterization  before  irrad- 
iation. As  discussed  in  a previous  section,  four-point 
resistivity  measurements  have  been  compared  with  Van  der  Pauw 
resistivity  measurements  for  calibration  purposes.  Mobility 
measurements  have  been  taken  from  one  slice  of  each  ingot 
and  the  values  are  very  near  average  values  given  in  Table  3. 

Note,  however,  that  the  mobility  in  the  Wacker's  sample  is 
some  131  too  low.  This  is  probably  due  to  the  effect  of 
mixed  conduction  in  this  very  high  resistivity  sample. 

Figure  22  shows  the  axial  resistivity  and  minority 
carrier  lifetime  of  ingot  VZ-072  after  cutting  it  into  slices 
of  approximately  1"  thickness.  Figures  23  through  28  show 
the  radial  uniformity  of  these  samples  as  obtained  from 
four-point  measurements.  Figure  29  shows  the  uniformity  of 
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TABLE  6.  Summary  of  Undoped  Sample  Characterization  Before  Irradiation. 


MURR  # 

Source 

P(U  pt) 
(fl-cm) 

P (V.d.P. 
(ft-cm) 

) UD 

(cm2/V-sec) 

Type 

T.P. 

Hall 

T 

( psec) 

0 

Monsanto  FZ 

7^9.0 

P 

- 

1 

Rockwell 

VZ-072-3 

2560. 

P 

110 

2 

VZ-072-3 

2953. 

P 

110 

3 

VZ-072-3 

3156. 

P 

110 

1* 

VZ-072-3 

2721. 

1*61*.  75 

P 

P 

110 

5 

Monsanto  FZ 

1690. 

1383.1* 

N 

N 

- 

6 

Monsanto  FZ 

1037.5 

lUo4 . 5 

N 

N 

- 

7 

Rockwell 

VZ-072-2 

2U50- 

2530  (ingot) 

P 

112 

8 

Rockwell 

VZ-083-1* 

7630. 

7729. 

52U.8 

P 

P 

_ 

9 

Rockwell 

VZ-083-5 

8189. 

8067. 

527.8 

P 

_ 

10 

Monsanto  CZ 

(after  irrad 
& anneal) 

3.222 

1267.0 

N 

N 

- 

11 

Monsanto  FZ 

1376.2 

11*73.5 

N 

N 

- 

ll* 

Topsil 

36-520-0-111 

P 

1200 

15 

Topsil 

36-520-0-111 

1*710- 

6U67 

(1"  ingot) 

P 

1200 

l6 

Topsil 

36-520-0-111 

U7U2. 

P 

1200 

17 

Wackers 

30661-6 

66512.  ? 

P 

530 

18 

T.I.  Lopex 
L1092 

U607 

P 

250 

19 

Topsil 

36-52O-O-III 

U9L2. 

1*60. 

P 

P 

1200 

20 

T.I.  Lopex 
L1092 

5195.9 

1*1*3.86 

P 

P 

250 

21 

Wacker 

30661-6 

19370. 

P 

530 

( trapping 

22 

UoaVoi" 

30661-6 

201*20. 

U27.0 

P 

P 

530 

A 
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TABLE  6. 

Summary  of  Unhoped  Sample  Characterization  Before 

Irradiation 

(Cont. ), 

MURR  t> 

Source 

P (U  pt) 
(Ji-cra) 

tfV.o.P.)  ,JD 

(Ji-co)  (cn2/V-se 

Type 

c)  T.P.  Hal 

1 

1 (yse-c 

23 

Similar  to  MURR 

10 

N 

- 

2k 

Wacker 

30661-6 

For  NAA 

P 

530 

25 

T.I.  Lopex 

L1092 

For  NAA 

P 

250 

26 

Topsil 

36-52O-O-III 

For  NAA 

P 

1200 

I. 

27 

Topsil 

36-52O-O-III 

52U8 

P 

1200 

28 

Monsanto  FZ 

15^0 

N 

- 

29 

Rockwell 

VZ-083-17 

12509. 

P 

30 

Rockwell 

VZ-083-18 

12796 

P 

31 

Rockwell 

VZ-083-19 

13326 

P 



32 

Rockwell 

VZ-083-20 

13687 

P 

33 

Rockwell 

VZ-083-21 

11*160 

P 

. 

Zk 

Rockwell 

VZ-083-22 

11*382 

P 

E i 

35 

T.I.  Lopex 

LI  09  2 

5095 

P 

36 

T.I.  Lopex 

L1092 

5587.6 

P 

- 

37 

Rockwell 

VZ-065-1 

For  NAA 

P 

38 

Rockwell 

VZ-065-2 

For  NAA 

P 

38-45 

Rockwell 

VZ-065  wafers 

(not  used) 

P 

46 

Hughes 

F0-1770602 

6183-6776 

( ingot ) 

P 

1200 

47 

Hughes 

F0 -1771101 

5610-8360 

( ingot ) 

P 

2000 
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TABLE  6.  Summary  of  Undoped  Sample  Characterization  Before  Irradiation  (Cont.). 


MURR  # 


Source 


P(l*  pt)  p ( V.d.P. ) yD  Type  T 

(fi-cm)  (ft-cm)  (cm2/V-sec)  T.P.  Hall  (yaec) 


4 8 

Hughes 

F0-1770602 

5785  (fast  neutron  irrad,  ) 

P 

49 

Hughes 

F0-1770602 

5675  ( isochronal  anneal ) 

P 

50-64 

Rockwell  Ga 
GZ-156-6  to 

doped  detector  material 

GZ-156-20  from  VZ  107  _ 

0.431-0.450,  [B]  = 5.42  x 10  c®*j 

P 

before  ( 

65 

Rockwell 

GZ-156-G 

For  NAA 

Fig.  52  p(x)  in  Ga  doped 

P 

66 

Rockwell 

VZ-083-50-1 

17510  p(  x ) for  [P]/[B]  = 10 

67-68 

T.I.  Lopex 
L1092 

Irradiated  and  sent  to  Ramdas, 
Purdue  for  optical  measurements 

69 

Rockwell 

VZ-083-50-2 

Fig.  53  p(x)  in  half  wafer  vs  whole  wafer 

70 

Rockwell 

VZ-083-50-3 

18970  p(  x ) for  [P]/[B]  = 1.12 

71 

Rockwell 

GZ-141-2 

p(x)  in  conventionally  compensated  sample 

72 

Rockwell 

GZ-141-3 

p(  x ) in  conventionally  compensated 

sample 

73 

Rockwell 

GZ-156-23 

For  NAA  in  Ga-doped  (not  completed) 

m 

Hughes 

F0-177062 

For  NAA  (not  completed) 

77 

Wackers 

102531-H<100>  p vs  1/T  (Ga  doped)  Fig.  32 

78 

Rockwell 

GZ-156-H 

p(x)  in  Ga  doped 

79 

Rockwell 

VZ-083-50-4 

19460  p( x ) for  [P]/[B]  = 2.04 

80 

Rockwell 

VZ-083-50-5 

19480  p(  x ) for  [P]/[B]  = 3.52 

81 

Rockwell 

VZ-083-50-6 

19380  p(  x ) for  [p]/[b]  = 4.45 

85 

Rockwell 

GZ-156-21 

p vs.  1/T  (Ga-doped)  Fig.  33 

86 

Hughes 

F0- 177060 2 

p(  x ) Fig.  30 

66 

TABLE  6.  Summary  of  Undoped  Sample  Characterization  Before  Irradiation  (Cont.), 


MURR  H 


To 

Check 

Uniformity 


Source 

p(4  Pt) 

(fl-cn) 

VZ- 072-1 

2270-2430 

VZ-072-2 

2450-2530 

VZ-072-3 

2520-2580 

VZ-072-4 

2580-2560 

VZ-072-5 

2550-2490 

(V.d.P.)  WD  Type 

(fi-cm)  (cm2/V-sec)  T.P.  Hall 


p(x)  Figs.  23-28 


T 

(psec) 


VZ-083,  while  Figure  30  shows  the  uniformity  of  a Hughes  sample. 

The  rise  in  resistivity  near  the  ends  of  these  wafers  is  due 
to  edge  effects  and  is  not  real.^ 

Table  7 summarizes  the  minority  carrier  lifetime  for 
the  various  ingots  which  have  been  measured. 


TABLE  7.  Summary  of  Minority  Carrier  Lifetimes. 


Source 

Ingot  No. 

T Measured 

T Specified  by  Manuf. 

Rockwell 

VZ-072 

100  psec 

100  psec 

Topsil 

36-52O-O-III 

1200  (trapping) 

1700 

T.I.  Lopex 

L 1092 

250 

500 

Wacker 

30661-6 

530 

3200 

Hughes 

FO-1770602 

1200 

900 

We  believe  that  the  lifetimes  in  this  resistivity  range  are 
probably  over  exaggerated  by  the  suppliers  because  of  over- 
modulation in  experimental  set  ups  designed  in  general  for 
lower  resistivity  material.  It  should  be  noted  that  our 
observed  lifetimes  for  the  Rockwell  material  are  in  general 
agreement  with  their  measurements.  The  lifetime  of  the  Topsil 
material  indicates  some  evidence  of  trapping  that  does  not 
saturate.  It  is,  therefore,  possible  that  the  real  lifetime 
is  shorter.  We  would  estimate,  however,  that  it  is  probably 
not  shorter  than  500  psec.  The  lifetime  of  the  Hughes  sample 
was  the  largest  found  and  was  free  of  trapping.  This  sample 
is  a [ill]  pull,  however. 

Table  8 summarizes  the  drift  mobilities  observed  for  the 
various  ingots  obtained  from  Hall  mobilities  using  the  relation- 
ship for  the  Hall  factor  shown  in  the  table. 

TABLE  8.  Summary  of  Drift  Mobilities. 


Source 

Ingot  No. 

Type 

yD  (cmVv-sec) 

Rockwell 

VZ-072 

P 

I465 

Rockwell 

VZ-083 

P 

525-528 

Topsil 

36-520-0-111 

P 

U60 

T.  I.  Lopex 

L1092 

P 

kkk 

Wacker 

30661-6 

P 

U27 

Monsanto  FZ 

- 

N 

1383-1^73 

Typical  (from 

Table  3 and  p.  3) 

P 

U9U 

N 

1396 

rp  = “h^d 

= 0.77 

rn  = 1-3 
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Center 


resistivity  of  VZ-072-5-bot 


30.  Radial  resistivity  of  a typical  Hughes 
wafer. 


We  conclude,  based  on  the  data  in  Tables  7 and  8 that  the 

materials  from  these  sources  are  of  comparable  crystal  quality. 

The  Rockwell  material  exhibits  inferior  lifetime  and  crystal 

perfection;  however,  the  material  was  pulled  in  a [100] 

direction,  which  is  more  difficult.  Apparently,  a [ill] 

Rockwell  pull  would  be  of  comparable  crystal  quality.  In 

addition,  the  Rockwell  material  has  the  advantage  of  extremely 

uniform  boron  distribution  as  evidenced  by  Figures  23-28. 

Carrier  concentration  from  Hall  effect  measurements  as  a 

function  of  reciprocal  temperature  have  been  measured  for 

several  representative  samples.  A Monsanto  Czochralski 

sample  which  had  been  heavily  transmutation  doped  was  first 

used  to  check  temperature  calibration  and  system  operation. 

These  data  are  shown  in  Fig.  31. 

The  solid  line  in  Fig.  31  is  a theoretical  fit  to  these 

data,  where  we  have  included  the  effects  of  the  first  excited 

2 1 

state  as  well  as  the  ground  state  of  phosphorus.  The  dotted 
curve  is  obtained  from  the  same  model  but  neglecting  the  first 
excited  state.  We  have  assumed  in  the  data  represented  in 
Fig.  31  that 

n = V 

The  theoretical  fits  are  obtained  from  the  relation 


(25) 
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2 2 

which  is  equivalent  to  the  more  usual  relationship 

n(n  + N.)  _ _ 

-TfD  - UA  -n  - * ■ 6Nce-ED/kT 

2 3/2 

where  N^,  = 2(2tt  m*kT/h  ) ' is  the  conduction  band  density 
of  states  and  8 = 1/2  is  a factor  representing  the  degeneracy 
of  the  phosphorus  ground  state  which  is  2.  We  have  assumed 
for  both  fits  to  the  data  that  Ep  = 0.044  eV  and  that 

7 7 - ? 

m*  = 1.08  m . To  take  the  statistical  effects  of  the 

e 

first  excited  state  into  account,  we  have  used  the  relation 
„ _ 2 ( 2tt  m*kT/h2)3/2  o‘E0/kT 

where  A = 0.010  eV  is  the  energy  difference  between  the  tenfold 

degenerate  first  excited  state  and  the  doubly  degenerate 

22 

ground  state  of  phosphorus  in  Eq.  (25).  Using  the  effect 
of  the  first  excited  state  gives  a slightly  better  fit 
than  neglecting  it. 

The  best  fit  to  the  data  determines  the  donor  and 
acceptor  concentration  as 

ND  - NA  = 1.2  x 1015  cm'3 

N.  * 1.8  x 1013  cm'3 
A 

and  a compensation  ratio  k = N^/Np  = 0.015.  The  acceptor 
concentration  is  a reasonable  number  for  the  boron  concen- 
tration in  Czochralski  Si. 
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The  deviation  from  the  fit  at  low  temperature  (i-  > 0.035) 
is  due  to  impurity  conduction.  This  effect  is  usually  seen  in 

? 7_  29 

compensated  material"  and  is  caused  in  Fig.  31  by  a carrier 

tunneling  from  an  occupied  donor  level  to  an  unoccupied  donor 

level.  Although  several  models  for  explaining  the  change  in 

activation  energy  at  low  temperatures  as  a function  of  com- 

27-29 

pensation  exist  for  resistivity  data,  no  adequate  model 

has  been  proposed  to  explain  the  Hall  coefficient  data.  ’ - 
It  is  not  clear  in  Figure  31  whether  the  carrier  concentration 
goes  through  a minimum  between  the  last  two  data  points  or 
whether  the  minimum  is  at  still  lower  temperature.  Since  the 
effect  of  impurity  conduction  occurs  at  higher  temperatures 
in  Si  than  in  Ge,  because  of  the  relative  depths  of  similar 
impurities,  it  is  doubtful  that  low  boron  concentration  in 
a gallium  doped  sample  can  be  determined  with  any  high  degree 
of  accuracy  because  of  interference  from  impurity  hopping 
conduction  induced  by  the  always  present  compensation. 

The  above  model  for  n-type  material  with  a single  donor 
level  and  acceptor  levels  far  removed  from  the  Fermi  level 
can  be  extended  to  the  case  of  Ga  doped  Si.  In  this  p-type 
material,  two  acceptor  levels  dominate  the  position  of  the 
Fermi  level  and  the  donor  levels  are  far  removed  from  it. 

To  describe  this  situation,  we  use  the  same  model  as  the 
Hughes  group.30 

Using  charge  conservation 

p * nd  ■ N1  *'CEp  * n2  £'CE2) 
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where  N£  is  the  concentration  of  gallium  atoms,  N2  the 
concentration  of  boron,  f'  (E)  the  Fermi  function  for 
hole  occupation  of  an  energy  level  at  energy  E above  the 


valence  band,  and  E£  and  E2  are  the  gallium  and  boron 


ground  state  ionization  energies  respectively.  Then,  using 


the  Boltzman  approximation 


p = N 


E£/kT 


V 


to  eliminate  the  Fermi  energy  from  the  above  equation,  we 
find  that 
p + N 


N, 


N, 


D rig  c*it) 


and  rearranging,  we  arrive  at  the  cubic  equation  for  the 
hole  concentration,  p, 

(p  ♦ ND) (p  + K:)(p  + K2)-N1K1(p  ♦ K2)-N2K2(p  + K£)  = 0 


where 


Ny  -E./kT 
— . e 1 

n ^ 


K2  = 


6 


V -E2/kT 
e 


Ny  = 2C2irm*kT/h2)3/2. 


We  have  used  the  following  parameters  to  obtain  fits  to  the 
experimental  data: 

E£  (gallium)  = 0.065  eV 

E2  (boron)  = 0.045  eV 

m*  = 0.59  m 

0 

B = 4 . 
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The  temperature  dependence  of  the  effective  mass  and  the 
acceptor  excited  states  have  been  neglected. 


A computer  code  to  solve  the  cubic  equation  for  each 
value  of  1/T  has  been  completed.  This  code  uses  the  "Half 
Interval  Search"  method  to  find  the  positive  real  root  of  the 
cubic  equation  for  p given  above. 

The  data  for  a Ga  doped  wafer  from  Wacker  are  shown  in 
Fig.  32  along  with  several  theoretical  fits  calculated  from 
the  above  model.  The  curves  are  not  particularly  sensitive 
to  and  NQ.  Our  best  estimates  for  these  concentrations  are 
Nx  = N(Ga)  = 1.78  x 1016em'3 

N2  = N (B)  = (1.25  + 0.25)  x 1013  cm' 3 

12  13-3 

1 x lO1^  = Nd  <<  1 x 10  cm  J . 

Also  included  are  curves  for  Ng/Ng  = 1.5  (curve  3). 

Ng/No  = 1.0  for  exact  compensation  (curve  4),  and  Ng/Ng  =3.0 
for  over  compensation  (curve  5).  It  is  clear  that  some  in- 
formation about  compensation  can  be  obtained  from  this  type 
of  data,  however,  the  sensitivity  to  changes  in  Ng  is  not 
very  good  for  cases  of  interest. 

Figure  33  is  a similar  carrier  concentration  plot  vs. 
reciprocal  temperature  for  a Ga  doped  Rockwell  sample, 
GZ-156-21.  These  data  have  not  been  analysed  in  detail  yet 
since  they  are  recently  taken. 

B.  Neutron  Activation  Analysis 

An  instrumental  neutron  activation  analysis  (INAA) 
procedure  for  measuring  trace  element  impurities  in  high 
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p (cm 


Figure  32.  Carrier  concentration  vs.  reciprocal 
temperature  for  a IVacker  Ga  doped 
float  zone  wafer  102531-11  (MURR-77) 


p (cm 


Figure  33.  Carrier  concentration  vs.  reciprocal 
temperature  for  a Rockwell  Ga  doped 
float  zone  wafer  GZ-156-21  (MURR-85) 
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purity  silicon  was  established.  This  INAA  procedure  serves 
as  a screening  technique  for  a wide  variety  of  unexpected 
impurity  elements  and  also  provides  quantitative  upper  limit 
data  for  a number  of  important  impurity  elements  if  they  are 
not  observed.  This  latter  group  of  important  elements  include 
(a)  elements  known  to  produce  striking  effects  on  the  physico- 
chemical properties  of  silicon,  (b)  elements  which  are  known 
to  be  common  impurities  in  technical  grade  silicon,  and  (c) 
elements  which  may  serve  as  indicators  for  contamination  of 
the  silicon  ingots  during  refining  and  handling. 

Short  Irradiation  Procedure.  A thick  wafer  weighing 
approximately  1 gram  was  etched  with  HF-HNOj  and  rinsed  with 
deionized  water  to  remove  any  surface  contamination.  The 
cleaned  sample  was  wrapped  in  a precleaned  polyethylene 

1 3 

bag,  and  irradiated  for  20  min.  at  a neutron  flux  of  8 x 10 
2 

n/cm  /sec  using  the  pneumatic  tube  transfer  facility. 

Following  irradiation  the  sample  was  washed  with  hot  8N  HNO^ 

and  rinsed  with  deionized  water.  The  gamma-ray  spectrum  of 

the  irradiated  sample  was  recorded  for  2000  seconds  after  a 

decay  period  of  60  minutes  using  a Ge(Li)  detector.  A three 

millimeter  Pb  absorber  was  placed  between  the  sample  and 

31 

detector  to  reduce  the  Si  beta  bremstrahlung.  Elemental 
standards  prepared  from  high  purity  elements  or  weighable 
compounds  were  irradiated  and  counted  in  the  same  manner  as 
the  samples. 

Long  Irradiation  Procedure.  Following  the  short 
irradiation  procedure,  the  same  samples  were  wrapped  in 


86 


aluminum  foil,  placed  in  a reflector  irradiation  canister, 


and  irradiated  for  40  hours  at  a neutron  flux  of  approximately 
14  2 

1 x 10  n/cm  /sec.  After  irradiation  the  samples  were  etched 
with  HF-HNOj  and  rinsed  with  deionized  water  to  remove 
surface  contamination  from  the  aluminum  foil  wrapping.  Each 
sample  was  counted  for  2000  sec  after  decay  periods  of  approx- 
imately 40  hours  and  5 days.  The  gamma-ray  spectra  from  0 
to  4 MeV  was  recorded  using  a Ge(Li)  detector.  Elemental 
standards  were  irradiated  and  counted  according  to  the  same 
schedule.  The  analytical  conditions  are  summarized  in  Table  9. 

Data  Reduction.  Existing  computer  codes  were  used  to 
reduce  the  gamma-ray  spectral  data  for  both  samples  and 
standards  to  peak  energies  and  peak  areas  for  all  observed 
gamma-rays.  Isotope  assignments  were  based  on  calculated 
energies  and,  where  necessary,  half-life  measurements. 

Elemental  concentrations  for  observed  impurities  were  cal- 
culated by  comparison  of  sample  and  standard  peak  areas  after 
proper  normalization  of  differences  in  irradiation  exposure 
and  decay  time,  if  necessary.  Minimum  detectable  peak  areas 
for  isotopes  not  observed  were  estimated  as  three  times  the 
statistical  deviation  of  the  spectral  base  line  in  the  peak 
region.  Detection  limits  for  the  elements  measured  are  given 
in  Table  9 for  a 1 gram  silicon  sample  using  a value  of  2.34 
gm/c.c.  for  the  density  of  silicon. 

Results . Silicon  ingots  obtained  from  four  different 
suppliers  were  analyzed.  Elemental  concentration  values 
are  given  in  Table  10  as  atoms  of  impurity  per  c.c.  of  silicon 
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TABLE  9.  Summary  of  Analytical  Conditions  and  Detection  Limits  for 
Neutron  Activation  Analysis  of  Silicon. 


Detection  Limit* 


Element 

Isotope 

E y(  keV ) 

Half-Life 

( atoms/cc ) 

20  min.  Irradiation  - 60  min. 

Decay  - 2000  sec 

Count 

In 

116m, 

In 

1293.1 

1 X 

io11 

Mn 

56Mn 

1812 

1.5 

x 1012 

40  hour  Irradiation  - 40  hour  Decay  - 2000  sec 

Count 

Mo 

"Mo 

140.6 

66.6  h 

4 x 

io12 

Cu 

64Cu 

511 

12.75  h 

7 x 

io11 

Ga 

72 

Ga 

834.1 

14.1  h 

3 x 

io11 

Au 

198Au 

411.8 

64.7  h 

2 x 

109 

W 

187w 

685.7 

23.9  h 

1 X 

io11 

As 

76  As 

559.2 

26.4  h 

1 X 

io11 

Sb 

122Sb 

564 

65.8  h 

6 x 

io10 

Yb 

175Yb 

396.1 

101  h 

8 x 

io10 

Na 

24Na 

1368.4 

15.0  h 

2 x 

io12 

Eu 

152mEu 

841.6 

9.3  h 

5 x 

io10 

40  hour  Irradiation  - 5 day  decay  - 2000  sec 

Count 

Cr 

51Cr 

320.0 

27.8  d 

9 x 

io12 

Th 

233Pa 

311.8 

27.0  d 

1.5 

, ,11 
x 10 

Sc 

46Sc 

889.4 

83.8  d 

1.5 

x 1011 

Hf 

18lHf 

482.0 

42.4  d 

3 x 

io11 

Co 

60Co 

1173.1 

5.25  y 

8 x 

io12 

Ag 

110m, 

Ag 

657.8 

253  d 

3 x 

io12 

Ta 

182Ta 

1121.2 

115  d 

1.5 

x 1012 

#3  a detection  limit  estimated  for  a 1 gram  sample  assuming  a density 
value  of  2.34  grams/cc  for  silicon. 
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TABLE  10.  Elemental  Concentrations  Measured  in  Four  Silicon  Ingots. 


Element 

Impurity  level  measured  in  atoms  of 

impurity/c.  c.* 

TI  Lopex  L1092 

Rockwell  VZ-072 

Topsil  36-520 

Wackers  30661-6 

In 

n.d. 

n.d. 

n.d. 

n.d 

Mn 

n.d. 

1? 

^2  x 10 

n.d. 

n.d. 

Mo 

12 

-V4  x 10 

n.d. 

n.d. 

n.d. 

Cu 

12 

2.9  x 10 

V7  x 1011 

12 

1.1  x 10 

12 

1.3  x 10 

Ga 

n.d. 

n.  d. 

n.d. 

n.d. 

Au 

4.7  x 109 

n.d. 

5.3  x 109 

^3  x 109 

As 

n.d. 

n.d. 

n.d. 

n.d. 

Sb 

n.d. 

n.d. 

n.d. 

n.d. 

Yb 

n.d. 

n.d. 

n.d. 

n.d. 

Na 

8.9  x 1013 

2.0  x 1013 

1.8  x 1013 

1.9  x 1013 

Eu 

n.d. 

n.d. 

n.d. 

n.d. 

Cr 

n.d. 

n.d. 

n.d. 

n.d. 

Th 

n.d. 

n.d. 

n.d. 

n.d. 

Sc 

n.d. 

n.d. 

n.d. 

n.  d. 

Hf 

n.d. 

n.d. 

n.d. 

n.d. 

Co 

n.d. 

n.d. 

n.d. 

n.d. 

Ag 

n.d. 

n.d. 

n.d. 

n.d. 

Ta 

n.d. 

n.d. 

n.d. 

n.d. 

W 

n.d. 

n.d. 

n.d. 

n.d. 

*n.d.  indicates  not  detected  concentrations  above  the  background. 

Detection 

limits  are  given  in  Table  9. 
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assuming  that  the  impurity  is  homogeneous  in  the  matrix.  The 

values  given  for  sodium  probably  represent  surface  contamina- 

27  24 

tion  derived  from  the  A1  (n,  a)  Na  reaction  occurring 
in  the  aluminum  foil  encapsulation. 

Figure  34  illustrates  a sawing  contamination  which  was 
detected  on  VZ-072.  This  sample  was  cut  on  a glass  cut-off 
tool  outside  our  laboratories.  It  is  apparent  that  gold 
surface  contamination  was  present  and  was  difficult  to  remove 
by  etching.  The  etching  procedure  apparently  tends  to 
selectively  etch  around  the  gold  spots  and  undercut  until 
the  gold  is  removed. 

C . Irradiation  and  Annealing  Experiments 

It  is  instructive  at  this  point  to  consider  the  types  and 
amounts  of  radiation  damage  which  result  from  a typical 
neutron  transmutation  irradiation.  We  will  consider  the 
following  sources  of  damage: 

• Fast  neutron  damage  for  a graphite  moderated  neutron 
energy  spectrum  (this  case  will  be  an  approximation  to 
our  sample  irradiation  position). 

• Fast  neutron  damage  from  a fission  neutron  energy 
spectrum  (this  case  will  approximate  in-pile 
irradiations  found  in  the  literature). 

• Gamma  damage  from  fission  gammas  from  the  reactor  pile. 

• Gamma  recoil  damage  (this  damage  occurs  as  a result  of 
conservation  of  momentum  when  a gamma  is  emitted  from 
an  excited  Si  nucleus  after  thermal  neutron  capture). 
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PERCENT  OF  SAMPLE  WGT.  REMOVED 


Figure  34.  Concentration  of  gold  vs.  percent  of 
sample  weight  removed. 
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• Beta  recoil  damage  (this  damage  occurs  as  a result  of 
conservation  of  momentum  when  the  beta  is  emitted  in 
the  deexcitation  of  "^Si  -+■  + fj~. 

We  will  neglect  damage  from  (n,  p) , (n,  a),  and  (y,  n) 
etc.,  reactions  because  of  their  low  cross  sections  for 
neutron  energies  of  interest.  Although  (n,  y)  reactions 
with  boron  have  large  cross  sections,  the  boron  concentration 
in  our  samples  is  negligible  compared  to  the  silicon  con- 
centration. Also,  we  estimate  that  the  damage  from  the  g‘ 
emitted  during  ^*P  production  will  produce  about  the  same 
number  of  displacements  as  the  recoil. 

Although  the  estimates  of  the  number  of  displacements 
due  to  each  of  these  processes  is  very  approximate,  it  is 
instructive  to  make  a rough  calculation  to  estimate  the 
magnitude  of  each  of  these  effects. 

The  number  of  displaced  atoms  per  unit  volume  per  sec  is 
calculated  from  the  equation 


^r—  = Nto<}>v  (2( 

where  Np  is  the  number  of  displacements  per  unit  volume, 

Nt  the  number  of  target  atoms  per  unit  volume,  <J>  is  the  flux 
of  damaging  particles  (#  particles/cm  /sec)  and  v is  the 
number  of  displacements  per  incident  damaging  particle.  We 
will  estimate  v from  the  Kinchin-Pease  model. ^ In  this 


model 


v(Er)  - ER/2Ed 


* 


for  ER  > 2E^  where  ER  is  the  average  kinetic  energy  of  the 
recoiling  lattice  atom  and  is  the  displacement  threshold 
in  silicon  which  we  take  as  about  12  eV.  The  maximum  recoil 
kinetic  energy  of  a target  atom  struck  by  a neutron  of  energy 


En  is  given  by 


4E. 


ER(max)  * — ^ 


n 


(28) 


when  A is  the  target  material  atomic  weight. 
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For  a fission  neutron  spectrum,  the  average  recoil 
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kinetic  energy  is  given  by' 


"R(fission  neutron) 


"R 


(max) 


(29) 


while  for  a graphite  moderated  spectrum 

1/2  ER(max) 


JR  (moderated  neutron) 


In 


ER(max) 

Lt; — J 


(30) 


We  shall  assume  that  the  average  neutron  energy  to  be  used 
in  equation  (28)  is  1.5  MeV.  Then  we  find  from  (28),  (29), 
and  (30)  that 


En (max)  = 4(1.5  MeV)  = 214  keV. 

R 2 i 


ER(fission  neutron) 


214  keV 
2 


107  keV 


and 


Er  (moderated  neutron) 


In 


/ 2 ( 2 1 4 keV) 


214  x 10J  eV 

— n~iv 


-] 


10.9  keV. 


We  find,  using  (27),  that  the  number  of  displacements 
per  neutron  is  much  higher  for  a fission  spectrum  than  for  a 
moderated  spectrum,  i.e., 


93 


v(fission)  = ^JT^) 

and 


4458  displacements/neutron  (31) 


10  9 keV 

v(moderated)  = 3 — = 455  displacements/neutron.  (32) 
Furthermore , in  our  bulk  pool  facility,  assuming  a Cd  ratio 
of  about  30:1,  the  fast  flux  is  about  (5  x 1011n/cm2/sec) /30 
* 1.67  x 10  n/cm  /sec.  For  an  in-core  irradiation,  the  Cd 
ratio  is  about  1:1  and  a typical  fast  flux  is,  therefore, 

5 x 10  n/cm  /sec.  Using  these  fluxes  and  assuming  that 
the  displacement  cross  section  is  approximately  equal  to  the 

7 •z  'J  A 

scattering  cross  section,  which  we  will  assume  is  3 x 10 
2 

cm  , we  find  using  (26)  and  (31) 

(fast  neutron  in-core)  = (5  x 1022) (3  x 10*24)(4458)  x 

3tT 

(5  x 102^)  = 3.35  x 1016  displacements/cm^/sec . (33) 

Using  (26)  and  (32)  we  find  also  that 


(fast  neutron  in-pool)  = (5  x 1022)(3  x 10'24)(455)  x 

(1.67  x 10*^)  = 1.14  x 1012  displacements/cm^/sec.  (34) 

The  gamma  flux  radiation  damage  is  also  quite  different 

since  in  core  the  photon  flux  is  about  3 x 106R/sec  while  in 

2 

our  bulk  pool  facility  it  is  about  3 x 10  R/sec.  The  cross 
section,  o(y)»  for  displacements  is  dominated  by  Compton 
electron  displacements  and  is  nearly  independent  of  photon 

<7  *r  7/1 

energy  in  silicon  for  typical  fission  gamma  energies.  ’ 

- 25 

We  will  assume  that  this  cross  section  is  about  1 x 10 
cm2  for  100  keV  < < 2 MeV  and  that  1 R/sec  = 2 x 10^ 

1 MeV  photons/cm2/sec. 35 


Then  the  number  of  displacements  is 
94 


(35) 


dND  (gamma  in  pile)  = (5  x 1022) (1  x 10_25)(3  x 106)x 

IF" 

(2  x 10^)  = 3 x 1013  displacements/cm3/sec 
and 

dND  (gamma  bulk  pool)  = (5  x 1022) (1  x 1C"25)(3  x 102)x 

IF" 

(2  x 10^)  = 3 x 10^  displacements/cm3/sec.  (36) 

The  various  silicon  isotopes  after  thermal  neutron 

absorption  relax  from  excited  states  by  prompt  gamma  emission. 

This  imparts  a recoil  to  the  silicon  isotope  in  order  to 

conserve  momentum.  Since 

Ev 

p = -I  = MV, 

E 2 

E = (1/2)MV2=  1/2  -X- 

K Me 

An  average  over  all  silicon  isotopes  and  cross  sections  yields 

37  38 

an  average  recoil  kinetic  energy  of  780  eV.  ' Therefore, 
using  (27),  we  find 

v = ZjJL = 32.5  displacements/thermal  neutron. 

2(12  eV) 

37  - 2 5 

The  average  cross  section  for  these  processes  is  1.3  x 10 
2 

cm  . The  flux  of  interest  is  the  thermal  neutron  flux,  there- 
fore 

dND  (y  - recoil  in  core)  = (5  x 1022) (1 . 3 x 10"25)(5  x 1013) x 

It- 

(32.5)  = 1.06  x 1013  displacements/cm3/sec  (37) 

and 

dNP  (Y-recoil  in  pool)  = (5  x 1022)  (1.3  x 1025)  (5  x 101 1 ) x 

It 

(32.5)  = 1.06  x 1011  displacements/cm3/sec . (38) 
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A similar  procedure  can  be  applied  to  beta  decay  recoil. 
Assuming  the  beta  energy  to  be  Eg  = l.S  MeV,  and  using 
conservation  of  momentum,  we  find 

P = /e62  - (mQC2)2  = MV 

c 


Er  = (1/ 2)MV2 


1/2 


C%c2)2 


T 


33.2  eV. 


Therefore,  using  (27), 

v = 2^2  eV)  = ^.76  displacements/neutron  absorbed. 

Since  this  occurs  for  only  31  of  the  total  silicon  atoms, 

- 2 5 2 

with  a cross  section  of  1.1  x 10  cm  , 

^ (B  recoil  in  core)  = (0.03) (5  x 1022)(1.1  x 10'25)  x 
(5  x 1013) (2 . 76)  = 2.3  x 1010  displacements/cm3/sec  (39) 


and 


^ (B  recoil  in  pool)  =»  (0.03)  (5  x 1022)(1.1  x 10"25)  x 

(5  X 10^*)  (2.76)  * 2.3  x 108  displacements/cm3/sec.  (40) 

* 

These  numbers  must  be  compared  with  the  rate  of  producing 
phosphorus  atoms  which  is 
dN 

^ (in  core)  = (0.03)(5  x 1022) (1 . 1 x 10‘25)(5  x 1013) 

■ 8.25  x 109  [P]/cm3/sec.  (41) 

and 

dN 

(in  pool)  = (0.03)  (5  x 1022)(1.1  x 10'25)(S  x 1011) 
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= 8.25  x 107 [P]/cm3/sec. 


(42) 


r 


We  can  now  compare  the  number  of  displacements  per 
phosphorus  atom  produced  by  each  of  these  mechanisms.  This 
information  is  summarized  in  Table  11. 

Although  these  calculations  are  quite  speculative, 
several  differences  between  in-core  (low  Cd  ratio)  and  in- 
pool (high  Cd  ratio)  irradiations  become  apparent.  It  is 
clear  that  the  fission  gamma  and  beta  recoil  damage  is  point- 
defect  like.  The  annealing  characteristics  of  this  damage 
should  be  similar  to  1 MeV  electron  irradiation  annealing. 

It  is  also  clear  that  the  gamma  recoil  damage  is  a significant 
fraction  of  the  radiation  damage  for  in-pool  irradiations 
but  is  relatively  less  significant  for  in-core  irradiations. 

Chukichev  and  Vavilov  have  estimated  experimentally  that 

the  number  of  gamma  recoil  displacements  is  about  equal  to 

the  number  of  fast  neutron  displacements  in  a heavy  water 

37 

moderated  reactor.  Our  ratio  of  fast  neutron  to  gamma 
recoil  displacements  is  of  the  order  of  10:1  for  an  in- 
pool position  with  a Cd  ratio  of  30:1.  Since  the  Cd  ratio 
in  a heavy  water  reactor  is  between  100:1  to  1000:1,  our 
calculations  are  in  reasonable  agreement  with  the  experimental 
data  in  Ref.  37. 

It  is  not  clear  whether  to  treat  the  gamma  recoil  radi- 
ation damage  as  point  defect  like  or  cluster  like.  The 
average  number  of  displacements  per  neutron  is  rather  large 
(32.5  displacements/thermal  neutron)  and  suggests  at  least 


TABLE  li.  Number  of  Displacements  per  Phosphorus  Produced. 


the  production  of  small  clusters.  The  actual  number  of 
displacements  per  thermal  neutron  absorbed,  however,  ranges 
from  about  10  or  20  up  to  80  to  160.  There  is,  therefore, 
a continuous  distribution  of  cluster  sizes  which  approach 
the  point  defect  distributions  at  the  lower  end.  The 
same  cluster  size  arguments  can  be  made  about  fast  neutron 
damage  from  a highly  moderated  energy  spectrum. 

We  have  considered  gamma  recoil  damage  both  as  cluster 
like  and  point  defect  like  in  Table  11,  to  calculate  the 
ratio  of  cluster  like  displacements  to  point  defect  like 
displacements  for  our  two  cases  of  interest.  It  seems  likely 
that  the  in-core  irradiations  deposit  3 to  80  more  dis- 
placements in  clusters  than  in  point  defects  as  compared  to 
in-pool  irradiations.  Such  a large  difference  in  the  defect 
spectrum  is  very  likely  to  produce  differences  in  annealing 
characteristics  for  the  two  cases.  We  will  discuss  some  of 
these  differences  in  the  following  discussion  on  annealing. 

Table  12  is  a summary  of  all  the  irradiations  performed 
thus  far,  the  sample  number,  the  fluence,  the  phosphorus 
concentration  added  by  irradiation  as  calculated  from  the 
cross  section,  and  initial  resistivity  assuming  a typical 
mobility,  and  the  purpose  of  the  irradiation.  In  this  table, 
4>c  is  the  fluence  required  to  take  the  particular  p-type 
sample  from  its  initial  resistivity  to  intrinsic  (or  exact 
compensation  of  the  residual  boron  by  the  phosphorus  added). 
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TABLE  12.  Summary  of  Irradiations  (undoped  material)  Cont. 
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MURR  29  " 6955  Approx  5000  ft-cm(N)  To  estatlish  limits  on  n-type  side 


MURR-70  RSAF  2565  3.975  x 10  estimated 

Rockwell  12-16 


An  inspection  of  this  table  shows  immediately  the  failure 
to  control  the  resistivity  for  partial  compensation  on  the 
p-type  side  of  intrinsic  (i.e.  $ < 4>c)  due  to  premature  type 
conversion. 

Figure  35  shows  this  premature  type  conversion  effect. 

This  figure  compares  an  isochronal  anneal  in  argon  of  a 

wafer  of  Rockwell  float  zone,  MURR  1,  with  the  annealing  in 

air  of  a piece  of  Monsanto  Czochralski,  and  vacuum  anneals 

38  39 

of  float  zone  and  Czochralski  by  Kharchenko  et  al.  * The 
solid  bar  is  the  expected  p-type  resistivity  for  the  Rockwell 
float  zone  after  the  final  anneal  on  the  basis  of  the  number 
of  donors  added  [0.042  ppb] . Instead  the  final  resistivity 
is  about  400  ft-cm  n-type  (the  type  conversion  probably  occurred 
between  600  and  800°C)  a gain  of  extra  donors  of  0.282  ppb 
(1.41  x 1013  atoms/cm3)  beyond  that  added  by  neutron  doping 
[0.042  ppb] . 

Also  shown  on  this  figure  is  the  isochronal  annealing  of 
a second  unirradiated  Rockwell  sample  from  the  same  wafer. 

It  should  be  noted  that  the  unirradiated  sample  shows  no 
adverse  effects  of  the  annealing  and  that  the  resistivity 
of  the  irradiated  sample  MURR  1 drops  below  that  of  the 
unirradiated  sample  at  300°C. 

The  heavily  irradiated  float  zone  sample  by  Kharchenko 
shows  no  type  conversion  peaks,  i.e.,  no  rise  in  resistivity 
toward  intrinsic  (230,000  n-cm) . This  is  confirmed  by  their 
Hall  measurements  directly  after  irradiation  but  before 

J 
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Resistivity 
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annealing.  It  should  be  noted  that  a drafting  error  has 

apparently  been  made  in  Figure  1 of  Ref.  39  since  the  final 

13-3 

values  of  n = 5 x 10  cm  and  p = 1 ft-cm  imply  an  electron 

2 

mobility  of  125,000  cm  /V-sec,  some  two  orders  of  magnitude 
too  high! 

Their  sample  was  apparently  irradiated  to  a fluence  of 
19  2 

about  2.68  x 10  n/cm  compared  to  the  Rockwell  sample 

16  2 

fluence  of  1.30  x 10  n/cm  . Therefore,  their  sample  ex- 
perienced a fluence  of  seme  2000  times  that  of  most  of  the 
samples  listed  in  Table  12.  Furthermore,  their  sample  was 
irradiated  in-core  while  ours  was  irradiated  in-pool  (see 
Table  11).  Because  of  the  fluence  differences  and  position 
differences,  it  is  estimated  that  the  number  of  displacements 
going  into  clusters  is  6 x 10^  to  6 x 10^  times  greater  than 
for  our  samples.  The  number  of  point  defects  is  about 
2000  times  greater. 

It  is  instructive  to  compare  the  radiation  damage  in 

the  float  zone  sample  by  Kharchenko,  et  al.  with  the  damage 

produced  by  a P+  ion  implanted  sample.  The  neutron  fluence 

for  a 1 ft-cm  n-type  final  resistivity  and  1000  fi-cm  p-type 

starting  resistivity  can  be  calculated  from  (7)  to  be 
19  2 

2.68  x 10  n/cm  . The  number  of  displacements  is  then,  using 

4 = <|>t,  and  equation  (26), 

ND  = (5  x 1022) (3  x 10'24) (2.68  x 1019) (4458) 

* 1.79  x 1022  displacements/cm^ . (43) 

22  3 

Since  silicon  contains  only  5 x 10  atoms/cm  , it  is 
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clear  that  the  radiation  damage  in  the  Kharchenko  sample 
is  approaching  the  amorphorus  condition. 

To  produce  an  amorphorus  layer  by  P+  ion  implantation, 
a dose  of  about  5 x 10^^  ions/cm^  at  40  keV  is  required. ^ 
The  number  of  displacements  per  ion  can  be  calculated  from 
(27)  to  yield 


E 

v = ' 7(TT"eV) 


40  keV  _ i_666  x io^  disp./ion. 


The  depth  of  a 40  keV  P+  implant  is  about  2 x 10  cm 
because  of  enhanced  diffusion  (see  Figure  2.27,  p.  56  in 
Ref.  40).  The  number  of  displacements  is 

= (5  x 1014  ions/cm^) (1 . 667  x 10^  displacements/ion) 
2 x 10  cm 


2 2 ^ 
= 4 x 10  displacements/cm  . 


(44) 


We  can  conclude  from  a comparison  of  (43)  and  (44)  that  the 
Kharchenko  float  zone  sample  is  nearly  amorphorus  as  a 
result  of  fast  neutron  damage. 

This  conclusion  has  interesting  consequences.  In  ion- 
implanted  silicon,  annealing  up  to  temperatures  of  about 

\ 

800°C  is  required  for  complete  electrical  activity  of  the 

implanted  ion  in  cases  where  an  amorphorus  layer  is  not 

formed.  The  recovery  of  the  electrical  activity  of  an 

amorphorus  layer  is  completed  by  about  600°C  and  is  closely 

related  to  the  formation  of  crystalline  films  of  Si  during 

4 3 

vacuum  evaporation  and  to  epitaxial  regrowth.  It  is 
reasonable  to  hypothesize  that  an  epitaxial  regrowth  of 
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cluster  damage  could  be  completed  by  600°C.  Further 
evidence  for  this  hyphothesis  will  become  apparent  later. 


Figure  35  also  shows,  for  comparison,  the  isochronal 
annealing  of  two  Czochralski  silicon  samples  (oxygen  con- 
centration  ^ 10“''cm  ~).  The  data  for  the  Kharchenko  sample 
is  taken  from  Ref.  39.  The  reverse  annealing  peaks  are 
evident  in  both  Czochralski  samples  and  are  apparently  re- 
lated to  the  high  oxygen  concentration  in  these  samples 
and  not  to  neutron  fluence  since  the  final  resistivities  differ. 


Evidence  for  this  point  of  view  is  shown  in  Figure  36. 

15  - 3 

These  float  zone  samples  (oxygen  concentration  'v.  10  cm  ) 
were  irradiated  to  a fluence  sufficient  to  produce  'vlOO  fi-cm 
n-type.  The  expected  final  resistivity  values  calculated  from 
the  starting  resistivity  are  shown  on  the  figure.  The  Topsil 
sample  came  to  within  1%  of  the  expected  donor  concentration 
while  the  Rockwell  sample  was  only  slightly  worse.  It  should 
be  noted  that  the  fluences  for  these  two  samples  are  similar 
to  the  fluences  for  the  Czochralski  samples  in  the  previous 
figure;  however,  the  annealing  peaks  and  dips  between  200 
and  500°C  do  not  appear  in  the  float  zone  material  at  this 
fluence.  Furthermore,  dislocation  density  does  not  appear 
to  have  a direct  correlation  with  the  lack  of  the  reverse 
annealing  peaks  shown  in  Figure  35  since  the  Rockwell  sample 
in  Figure  36  is  highly  dislocated  while  the  Topsil  sample  is 
not.  We  conclude  then,  on  the  basis  of  circumstantial  evidence, 
that  the  reverse  annealing  peaks  are  a result  of  the  concen- 
tration of  oxygen  relative  to  the  level  of  phosphorus. 
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Figure  36.  Isochronal  annealing  of  two  float  zone 
samples  irradiated  to  produce  ^100  fi-cm 
n- type . 
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Kharchenko  et  al.  also  state  that  their  Czochralski 

39 

samples  were  p-type  up  to  600°C  and  n-type  thereafter. 

(Their  Czochralski  samples  are  apparently  irradiated  in-pool.) 
We  conclude  from  our  data  in  Figure  36  that  electrical  activ- 
ity of  the  phosphorus  begins  at  about  600°C  in  agreement  with 
their  conclusions.  The  large  reverse  annealing  peak  from  600- 
800°C  is  primarily  associated  with  the  electrical  activation 
of  phosphorus.  The  Rockwell  sample  resistivity  is  intrinsic 
at  700°C  while  the  intrinsic  resistivity  peak  in  the  Topsil 
sample  has  apparently  been  partially  missed.  It  can  further 
be  concluded  from  the  data  in  Figure  36  that  nearly  all  the 
electrically  active  radiation  damage  has  annealed  at  temper- 
atures below  600°C.  This  may  indicate  the  completion  of  the 
annealing  of  clusters  in  agreement  with  the  hypothesis  pre- 
sented earlier.  This  contention  can  be  justified  as  follows. 

We  can  calculate  the  number  of  electrically  active  donors 
as  a function  of  isochronal  annealing  temperature  assuming  the 
mobility  has  completely  recovered  before  6>'0°C  and  the  reverse 
annealing  peak  between  600  and  800°C  is  the  result  of  a type 
conversion  from  p-type  to  n-type  due  to  the  electrical 
activation  of  phosphorus.  (The  thermal  probe  type  as  a 
function  of  annealing  temperature  is  shown  on  the  figure.) 

On  the  basis  of  this  assumption  we  can  calculate  the  fraction 
of  phosphorus  which  is  not  electrically  active  (N/Nq)  as  a 
function  of  annealing  temperature.  We  define  that  this 
fraction  is  unity  at  the  resistivity  minimum  for  normal- 
ization purposes.  Then  if  the  electrical  activation  is 
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governed  by  first  order  annealing  kinetics, 
d(N/No)  = -K(T)(N/N0) 

si — 

where  K(T)  = v Qe  is  the  temperature  dependent  rate  at 

which  the  fraction  of  phosphorus  not  electrically  active 
decreases  at  the  various  annealing  temperatures.  From  the 
above 

In  (Nq/N)  = (voAx)e'E/kT 

where  Ax  is  the  annealing  time  at  temperature  T.  For  equal 
annealing  times  at  various  temperatures,  a semilog  plot  of 
In  Nq/N  vs.  1/T  yields  the  activation  energy  for  the  process 
while  the  intercept  yields  the  defect  jump  frequency. 

The  result  of  this  first  order  analysis  is  shown  in  Figure 
37.  It  is  clear  that  the  first  order  annealing  assumption  is 
reasonable  and  that  the  first  order  annealing  occurs  through- 
out the  type  conversion  peak  in  Figure  36.  To  have  obtained 
this  result,  the  two  assumptions  mentioned  previously  must  be 
valid,  i.e.,  a lack  of  electrically  active  damage  at  tempera- 
tures above  600°  and  the  complete  recovery  of  carrier  mobility 
before  600°C. 

The  first  order  annealing  analysis  suggests  that  the 

activation  energy  for  the  phosphorus  to  become  electrically 

active  is  0.96  +_  0.05  eV  while  the  jump  frequency  is  about 
2 _ ^ 

2 x 10  sec  . This  activation  energy  is  remarkably  close  to 
the  values  given  in  the  literature  for  the  reorientation  of 
the  vacancy  phosphorus  E-center  (E  = 0.93  + 0.05  eV, 
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Figure  37.  Semilog  plot  of  the  In  of  the  reciprocal 
of  the  fraction  of  phosphorus  not 
electrically  active  vs.  reciprocal 
annealing  temperature  for  the  Topsil 
sample  shown  in  Figure  36. 


■it  .144 

v = 1.6  x 10  sec  ) and  for  the  electrical  disappearance 
o 

of  a deep  acceptor  (Ec-0.4)eV  thought  to  be  the  E-center 

(E  = 0.94  eV,  v = 108sec'1).45 
o 

Watkins  and  Corbett  suggest  that  the  lower  jump  frequency 
in  Ref.  45  indicates  the  diffusion  of  the  E-center  through 
the  lattice  over  great  distances  (10  /10  = 10  jumps  before 

dissociation  at  a vacancy  sink).  They  point  out  that  an  E- 
center  reorientation  followed  by  a vacancy  phosphorus  inter- 
change is  equivalent  to  an  E-center  jump  through  the  lattice. 
If  we  make  this  same  assumption,  then  we  are  observing 

1013/2  x 10^  = 5 x 10^°  jumps  before  annihilation!  For  a 

10  ° 
random  walk  of  5 x 10  jumps  of  a distance  of  * 2A,  the 

mean  distance  traveled  would  be 

Lq  = \Js  x 1010  (2  x 10-8cm)  - 4.5  x 10'3cm. 

This  is  of  the  order  of  the  smallest  sample  dimension  and 
would  suggest  vacancy  anihilation  at  the  sample  surface.  But 
this  can  not  be  the  correct  mechanism  for  E-center  annihila- 
tion because  it  would  deposit  all  the  phosphorus  at  the 
sample  surface.  We  know,  however,  that  phosphorus  is 
uniformly  distributed  in  transmutation  doped  silicon.  We 
must  conclude,  then,  that  free  migration  of  E-centers  in 
NTD-Si  does  not  occur  for  some  reason.  One  possibility  is 
that  the  E-center  becomes  pinned  in  some  way  to  its  lattice 
position  by  the  clustering  with  other  defects.  Since  there 
are  approximately  three  displacements  per  beta  emission  when 
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the  phosphorus  atom  is  created  by  transmutation,  it  is  possible 
that  a di-interstitial  or  di-vacancy,  E-center  complex  is 
formed  which  remains  stable  up  to  the  higher  annealing 


temperatures  observed  here. 

It  should  be  noted  in  Fig.  36  that  at  the  minimum  re- 
sistivity near  600°C  in  the  Topsil  sample,  the  hole  concen- 
tration is  about  6 x 10^4cm  . Since  the  initial  hole  con- 

centration was  about  2 or  3 x 1011cm  3,  we  can  safely  infer 
that  a very  large  acceptor  concentration  has  been  created  as 

a result  of  irradiation  and  annealing  to  550°C  (an  alternate 

14  - 3 

interpretation  would  be  the  destruction  of  6 x 10  cm 
donors  but  this  would  imply  a compensation  ratio  of  0.9996 
before  irradiation  which  is  impossible  because  of  the  initial 
resistivity) . 

It  will  be  shown  shortly  that  similar  acceptor  defect 
concentrations  have  been  observed  in  all  float  zone  studied 
to  date  and  that  this  defect  concentration  is  nearly  independ- 
ent of  irradiation  fluence  over  a fluence  range  of  50:1.  We 
conclude,  therefore,  that  the  magnitude  of  this  concentration 
is  dependent  on  the  properties  of  the  silicon  before  irra- 
diation. The  most  likely  impurity  to  be  found  in  highly 
zone  refined  silicon  in  a concentration  range  larger  than 
3 x 1013cm  3 (Rockwell,  Wacker)  to  5 x 1014cm'3  (Topsil  and 
T.I.)  is  oxygen.  We,  therefore,  suggest  that  these  acceptor 
defects  are  single  or  multiple  defect-oxygen  complexes,  i.e. 
the  oxygen  impurity  concentration  is  decorated  as  a result 
of  the  process  of  radiation  damage  and  annealing  to  550°C. 
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Figure  38  shows  the  recovery  of  resistivity  and  minority 
carrier  lifetime  on  a Rockwell  sample  1"  dia  x 1"  long.  The 
fluence  for  this  irradiation  was  50%  of  that  necessary  to 
produce  type  conversion.  The  annealing  shows  several  unusual 
features  which  are  found  only  in  the  Rockwell  samples  and 
which  could  possibly  be  related  to  the  dislocation  density. 

The  resistivity  between  200°C  and  800°C  is  multivalued. 

Shortly  after  removal  from  the  annealing  furnace  and  quenching 
to  room  temperature  in  trichloroethylene,  the  resistivity  is 
a minimum  but  rises  to  the  second  higher  value  while  storing 
the  sample  in  the  dark.  Upon  illumination  from  a tungsten 
light,  the  resistivity  returns  to  its  t = o value.  The 
resistivity  can  be  cycled  any  number  of  times.  The  reciprocal 
of  the  minority  carrier  lifetime  follows  the  resistivity  in 
general.  Very  heavy  trapping  is  observed  in  the  same  region 
as  the  largest  change  in  resistivity.  This  sample  type- 
converts  at  about  800°C.  It  should  also  be  noted  that  a 
plateau  in  the  recovery  around  300°C  is  similar  to  the 
Czochralski  sample  annealing  seen  previously. 

In  an  attempt  to  determine  if  sample  contamination  is 

responsible  for  the  unexpected  type  conversion,  a series  of 

isochronal  anneals  were  performed  in  argon  and  a vacuum 
_ 7 

of  5 x 10  torr.  The  fluences  for  these  irradiations  were 
in  the  range  of  f = 0.36  $>c  to  4>  = 0.56  <J>c  where  <t>c  is  the 
fluence  required  for  exact  compensation.  The  results  are 
shown  in  Figures  39  through  45.  A comparison  of  a pair  of 
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Figure  38.  Isochronal  annealing  of  resistivity  and 
minority  carrier  lifetime  of  Rockwell 
float  zone  (argon  anneal)  bulk  sample. 
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Figure  39.  Argon  isochronal  anneal  of  Rockwell 
float  zone  wafer. 
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vacuum  and  argon  anneals  are  shown  in  Figures  46  and  47. 

It  can  be  seen  that  all  of  these  samples  type  convert  at 
temperatures  between  600  and  900°C.  The  expected  resistivities 
calculated  from  the  concentration  of  phosphorus  added  is  also 
shown  on  this  series  of  figures. 

Figure  48  shows  an  argon  isochronal  anneal  for  a Hughes 
wafer  in  which  the  fluence  is  slightly  larger  than  the  previous 
irradiations.  For  this  sample  the  fluence  was  just  under  the 
critical  fluence,  $c.  Again,  we  see  that  type  conversion 
occurs  due  to  unexpected  donors. 

Note  that  the  minimum  occurs  at  the  same  temperature 
in  all  of  these  anneals  and  that  this  temperature  minimum 
is  identical  in  the  heavily  irradiated  samples  shown  in 
Fig.  36.  This  minimum  resistivity  varies  between  200  to 
2000  fi-cm  in  all  these  samples  justifying  our  estimates  of 
the  concentration  of  acceptor  defects  discussed  previously. 

An  unusual  feature  seen  in  most  of  these  anneals  are  the 
Czochralski  like  reverse  anneals  at  temperatures  below  600°C. 
These  reverse  anneals  are  most  pronounced  in  those  samples 
which  show  the  lowest  minimums  at  550  to  600°C.  This  is 
evident  in  Figures  40  and  41.  They  are  much  less  pronounced 
in  the  Rockwell  and  Wacker  silicon  which  we  suspect  indicates 
about  an  order  of  magnitude  less  oxygen  in  these  samples. 

The  final  resistivities  are  completely  uncontrollable 
for  these  light  irradiations.  Although  it  is  tempting  to 
suspect  sample  contamination  during  the  anneal,  we  have  been 
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Figure  46.  Comparison  of  vacuum  and  argon  isochronal 
anneals  for  Topsil  wafers. 
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Figure  48.  Argon  isochronal  anneal  of  Hughes  wafer 
where  <t>  « <t>c. 
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unable  to  detect  any  contamination  by  neutron  activation  ana- 
lysis. Furthermore,  Figures  46  and  47  suggest  there  is  very 
little  difference  between  argon  and  vacuum  annealing.  Figure 

40  shows  a type  conversion  at  600°C.  It  is  difficult  to 

12  3 

imagine  diffusing  3 x 10  atoms/cm  into  this  sample  from 
surface  contamination  at  600°C.  We  have  seen  evidence  of 
surface  contamination  for  some  anneals  above  950°C,  however, 
impurity  diffusion  is  much  more  likely  at  this  temperature. 

Figure  49  shows  isochronal  annealing  in  argon  for  an 
irradiation  of  $ = 1.02  $c,  i.e.  2%  overcompensation.  The 
low  resistivities  at  temperatures  below  300°C  were  the  result 
of  remounting  this  sample  with  wax  on  a hot  plate  after 
irradiation  for  relapping  after  etching.  It  is  seen  that  the 
final  resistivity  is  very  close  to  the  expected  value  except 
for  the  last  annealing  point  which  may  be  surface  contamination. 
Because  of  slight  overcompensation,  type  conversion  is  ex- 
pected in  this  sample.  The  etching  and  lapping  may  have  re- 
moved a sodium  surface  contamination;  however,  the  same 
procedure  on  other  samples  has  made  no  improvement  in  our 
ability  to  prevent  type  conversion. 

Fernaps  the  additional  phosphorus  concentration,  the 
extra  radiation  damage  or  the  final  charge  type  has  suppressed 
the  unknown  donor  which  causes  type  conversion  in  some  way. 

In  the  following  table,  we  summarize  some  of  the 
prominent  annealing  features  common  to  all  samples  investigated 
to  date.  The  parameters  listed  below  are  the  two  "Czochralski- 
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Figure  49.  Argon  isochronal  anneal  of  Topsil  wafer 
irradiated  to  4 * 1.02$c. 
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like"  peaks,  the  minimum  resistivity  temperature  and  the 
ratio  of  the  minimum  resistivity  to  the  before  irradiation 
resistivity . 

In  order  to  study  the  radiation  damage  effects  of  the 
fast  neutrons  without  the  radiation  effects  of  thermal  flux 
and  transmutation  doping,  two  samples  were  placed  in  boron 
shielded  cans  for  irradiation.  Because  of  boron's  extremely 
large  capture  cross  section  for  thermal  neutrons,  the  thermal 
flux  seen  by  the  Si  is  greatly  reduced  but  the  fast  neutron 
flux  penetrates.  These  two  samples  were  given  irradiations 
of  10,000  and  30,000  integrator  counts  respectively  in  the 
bulk  pool  variable  flux  facility  (RSAF) . The  fast  flux 
irradiations  correspond  to  fast  neutron  fluences  which  would 
accompany  thermal  fluences  of  approximately  <p  and  3*. 

The  fast  to  thermal  ratio  could  be  roughly  estimated  by  in- 
cluding an  A1 : Co  flux  wire  in  each  can  at  the  Si  position. 

The  thermal  to  fast  ratio  for  RSAF  was,  therefore,  determined 
to  be  between  20:1  to  40:1.  This  agrees  with  our  previous 
estimate  of  30:1. 

The  argon  isochronal  annealing  of  these  two  samples  is 
shown  in  Fig.  50.  It  is  apparent  that  the  annealing  of  these 
samples  is  quite  similar  to  the  thermal  neutron  irradiation 
anneals  presented  earlier.  We  must  conclude,  therefore,  that 
most  of  the  electrically  active  radiation  damage  is  caused  by 
fast  neutron  damage. 


130 


RESISTIVITY  (.ft -CM) 


Figure  50.  Argon  isochronal  anneals  of  fast  neutron 
damage.  Thermal  neutron  flux  was 
suppressed  by  using  boron  irradiation  can. 


Except  for  a narrow  temperature  range  around  300°C,  the 
annealing  of  both  samples  up  to  550°C  is  independent  of 
fluence.  In  the  region  from  550°C  up  to  the  highest  anneal 
temperature,  changes  in  fluence  produce  differences  in  the 
two  anneals  and  is,  therefore,  fluence  dependent. 

It  should  be  noted  that  both  samples  recovered  closely 
to  their  original  values  in  contrast  to  those  samples  which 
have  had  phosphorus  added  by  transmutation.  There  is  no  ten- 
dency in  either  of  these  samples  to  type  convert!  The  resis- 
tivity in  the  region  of  600  to  800°C  corresponds  to  a donor 
concentration  which  is  at  least  two  orders  of  magnitude  lower 
than  in  those  samples  containing  transmutation  produced 
phosphorus  for  light  irradiations.  The  phosphorus  is,  there- 
fore , dominating  the  compensation  process  in  this  range. 
Also,  it  should  be  noted  that  two  small  recovery  stages  in 
the  range  of  600  to  700°C  are  just  visible  in  these  two 
samples.  These  annealing  stages  were  completely  hidden  in 
1 the  transmutation  doped  samples.  (A  hint  of  one  of  these 

stages  is  seen  in  Figs.  41,  43,  47,  and  48.) 

From  all  of  the  above  experiments,  it  is  apparent  that 
a single  temperature  annealing  procedure  can  be  devised. 

This  procedure  is  as  follows: 

1.  The  annealing  boat  consists  of  a quartz  boat  with  a 

3/4"  dia  x 2"  long  silicon  support  ingot  sliced  in  half 
lengthwise.  This  assembly  remains  inside  a closed 
spectrasil  annealing  tube  and  is  prepared  by  preheating 
the  spectrasil  annealing  tube  and  boat  at  800°C  with 
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flowing  argon  through  the  tube  for  two  or  more  hours. 

(See  Fig.  51.) 

2.  A 2"  float  zone  wafer  is  cleaved  diametrically,  cleaned 
with  soap  solution  using  swabs,  D.I.  water  rinsed,  blown 
dry  with  argon  and  swabbed  with  trichloroethylene.  These 
semicircular  liners  are  used  to  protect  the  test  wafer 
from  contamination. 

3.  The  prepared  spectrasil  annealing  tube  is  cooled  and 
opened;  the  boat  is  moved  to  the  mouth  of  the  tube 
but  NOT  removed. 

4.  Float  zone  semicircles  (liners)  are  etched  in  CP-4 
(fresh)  solution  for  30  seconds  in  a dedicated  container 
and  rinsed  in  D.I.  water. 

5.  Tweezers  are  cleaned  in  Caro's  acid  (1:1  -I^C^  iH-^SO^) 
(fresh),  D.I.  water  rinsed,  and  argon  blown  dry. 

6.  Liners  are  removed  from  rinse,  argon  blown  dry  (holding 
with  tweezers),  and  placed  vertically  in  boat.  (This 
operation  is  done  close  to  the  mouth  of  the  tube). 

7.  The  sample  is  swabbed  with  soap  solution,  D.I.  water 
rinsed,  argon  blown  dry,  and  swabbed  with  tri- 
chloroethylene . 

8.  Tweezers  are  recleaned  as  in  #5. 

9.  The  sample  is  etched  in  CP-4  for  30  seconds  (dedicated 
container)  D.I.  water  rinsed,  argon  blown  dry  (while 
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Figure  51.  Annealing  furnace  tube,  sample  boat, 
and  silicon  spacers. 
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holding  with  tweezers)  and  placed  between  liners. 


10.  The  sample  boat  is  slid  into  the  annealing  position 
in  the  spectrasil  annealing  tube  such  that  the  liners 
are  flush  against  the  sides  of  the  sample  and  the 
spectrasil  tube  is  sealed  and  flushed  with  argon 

at  the  rate  of  about  30  cc/min. 

11.  The  spectrasil  annealing  tube  is  inserted  into  the 
furnace  at  800°C,  remains  for  one  hour,  then  power  is 
disconnected.  The  sample  cools  in  the  furnace  in 
flowing  argon  (15  cc/min.)  until  the  temperature  is 
less  than  100°C. 

12.  The  sample  is  removed,  lapped  with  1000  grit  using 
plastic  gloves  for  handling  and  D.I.  water  rinsed. 

The  reason  for  the  final  lap  after  annealing  is  to 
remove  shallow  impurity  concentrations  which  have  diffused 
into  the  surface  as  a result  of  handling.  The  importance 
of  this  is  shown  in  Figure  52.  The  annealed  sample  data 
points  represented  by  circles  show  the  effect  of 
accidental  touching  of  the  sample  before  anneal.  The 
lapping  removes  these  impurities  and  dramatically  improves 
resisitivity  fluctuations.  The  final  average  value 
near  800  fi-cm  is  very  close  to  the  expected  resistivity 
value  for  this  sample. 
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oistance  from  probe  center  to  wafer 


In  order  to  evaluate  the  success  of  our  single  anneal 
temperature,  it  is  useful  to  extend  the  models  presented  in 
Section  II  to  answer  the  following.  Given  the  initial 
resistivity  pQ  (p-type)  and  the  final  resistivity  p^  (n- 
type) , what  is  the  concentration  of  phosphorus  added  if  the 
number  of  intrinsic  carriers  is  not  negligible  and  a mixed 
conductivity  condition  exists? 

The  initial  carrier  concentration,  pQ,  is  given  by 

Po  - 1/2 [N  - Nd  * V (Na  - Nd  )2  . 4n  2'  ] 

O 0 

where  n^  is  the  number  of  intrinsic  carriers.  Solving  for 


(45) 


Likewise, 

• na  * y <NDf  - na>2  * 4niz  i 


and 


NDf  - NA  = 


nf  - ni 


n. 


(46) 


Adding  (45)  and  (46),  we  find  that  the  concentration  of 


phosphorus  added,  Np  = Nn  - Nn  , is 

y uf  uq 

2 2 2 2 

ND  . P-°—  - ' "i 

P P_  nr 


(47) 


This  equation  allows  a calculation  of  the  concentration  of 
phosphorus  added  by  transmutation  from  the  initial  number  of 
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holes  (p-type)  and  final  number  of  electrons  (n-type) . But 
these  concentrations  are  related  to  the  measured  resis- 
tivity and  b = yn/yp  by  Eq.  (13) 


= < 


7 7-  , P-type 

eyp  (bni  + p ) 


n 


n-type . 

Solving  these  relations  for  n or  p as  a function  of  the 


5 2 9 

^ eyp  (bn  + n/  ) 


(48) 


(49) 


Thus,  given  pQ  and  Pp,  we  can  calculate  pQ  and  n^  from  (48) 
and  (49),  and  by  substitution  into  (47),  we  obtain  the  con- 
centration of  phosphorus  added  in  the  presence  of  mixed 
conduction  for  the  case  where  n^  is  not  negligible.  It 
should  be  noted  that  (47)  reduces  to  Np  = pQ  + n^  if  a<<pQ, 
nf.  Likewise,  under  similar  conditions  (48)  and  (49)  become 


n. 


pfeyn 


* Po  " P„ey 


(50) 


For  p >>n.  but  nr  ^ n.  we  have  a condition  which  exists  for 
o l r l 

slightly  overcompensated  samples.  For  these  conditions 
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where 


p = ^ and  nr  is  obtained  from  (48)  . 

0 poeyp  f 

The  results  of  the  single  8 00°C  anneal  temperature 
procedure  are  show  in  Table  14.  The  success  of  achieving 
the  correct  resistivity  for  overcompensated  material  has 
been  as  spectacular  as  the  failures  for  undercompensation. 
In  Table  14,  the  initial  and  final  resistivities,  the 
concentration  of  phosphorus  added  as  expected  from  the 
integrator  calibration,  Np(E),  the  concentration  of 
phosphorus  added  as  calculated  from  the  resistivities, 
Np^  = K4>,  and  the  percentage  error  between  expected  and 
observed  phosphorus  concentration  is  shown.  Note  that 
a resistivity  of  over  100,030  S2-cm  has  been  achieved  in 
MURR-70.  This  represents  an  uncompensated  carrier  con- 
centration of  only  3.77  x lO^cm  about  three  times  the 
intrinsic  carrier  concentration.  The  conventional 
compensation  ratio  N^/Np  is  about  0.95.  Because  of  the 
success  of  these  last  experiments,  it  was  felt  that  the 
Ga  doped  Si  could  be  successfully  compensated.  The 
remainder  of  the  compensation  experiments  were,  therefore, 


j 

performed  on  samples  taken  from  Rockwell  ingot  GI-156, 
a Ga  doped  ingot. 


1 PPD  = 4.996  x 10 1 3cm" 3 
♦Argon  anneal  for  1 hr.  at  800°C  ♦ slow  cool. 

♦♦Calculated  approximately  from  initial  resistivity  before  1st  irradiation. 


determine  the  boron  compensation.  After  multiple  zone  passes 
to  reduce  impurities  other  than  boron  to  negligible  concent- 
rations but  before  the  addition  of  Ga,  a fourpoint  probe 
measurement  of  the  starting  ingot,  VZ-107,  indicated  that 
p = 2.4  x 10^  fi-cm  and  that  p = 0.108  ppb  = 5.42  x lO^cm 

This  value,  obtained  from  Rockwell,  was  used  to  determine  the 
critical  fluence,  <J>  . Ga  was  then  added  to  the  ingot. 

Although  this  procedure  should  have  worked  well,  several 
problems  arose.  The  current  and  voltage  leads  on  the  Rockwell 
apparatus  were  hooked  to  the  wrong  instruments  leading  to  an 
estimated  error  of  20%.  Furthermore,  a measurement  of  p(x) 
vs.  x on  undoped  wafers  from  this  ingot  (GZ-156  -C,D,G,  and 
H)  indicated  surface  contamination  from  unknown  sources 
This  is  shown  quite  clearly  in  Fig.  53  for  GZ-156-G.  The 
resistivity  near  the  edge  drops  to  about  1 fi-cm.  Since  this 
sample  was  only  a half  wafer  (its  shape  is  shown  in  the  upper 
right  of  Fig.  53).  the  edge  effect  on  the  four  point  probe 
technique  was  claimed  to  be  the  source  of  Lhe  difficulty  by 
the  Rockwell  group.  We,  therefore,  performed  a similar 
experiment  on  a half  wafer  from  VZ-083  both  before  and  after 
cleaving  into  a half  wafer  shape.  The  results  of  this  ex- 
periment are  shown  in  Fig.  54.  As  logic  and  theory  would 
dictate,  the  apparent  resistivity  for  the  half  wafer  from 
VZ-083  rises  at  the  edges.  We  can  only  conclude,  then,  that 
surface  contamination  exists  on  wafer  GZ-156-G.  We  have 
returned  this  wafer  to  Rockwell,  at  their  request,  for  further 
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DISTANCE  FROM  EDGE  OF  WAFER  TO  CENTER  OF  PROBE  (inj 


analysis.  A second  half  wafer  with  similar  p(x)  vs.  x 
characteristics  has  been  retained  by  MURR  for  activation 
analysis . 

Since  it  was  not  clear  which  value  of  overcompensation 
of  the  boron  concentration  would  produce  the  best  D*,  com- 
pensation ranges  for  4>/<t>c  from  about  1:1  up  to  20:1  were 
performed.  This  range  of  doping  ratios  was  requested  by 
Rockwell.  These  irradiations  are  summarized  in  Table  15. 

All  of  these  samples  were  slightly  underirradiated 
(-3  to  -10%)  due  to  a lack  of  calibration  factor  for  the 
high  number  of  integrator  counts  encountered  for  these 
samples.  This  is  of  little  consequence  because  of  the  204 
error  in  initial  resistivity  determination  and  the  nonuniform- 
ity of  resistivity  for  GZ-156-G.  One  set  of  samples  was  irra- 
diated, cleaned,  and  annealed  at  MURR  while  an  identical  set 
was  irradiated  but  not  annealed  or  cleaned.  Both  sets  were 
sent  to  Rockwell.  Three  samples  were  retained  by  MURR  for 
isochronal  annealing  experiments  (MURR- 52,  MURR- 53,  and 
MURR-56)  corresponding  to  d>/<|>  values  from  1 to  3.  All 
other  samples  from  this  set  were  returned  to  Rockwell  for 
detector  fabrication. 

Isochronal  anneals  of  the  Ga  doped  samples  (MURR-S2, 
MURR-53,  and  MURR-56)  were  performed  to  evaluate  any  new 
difficulties  which  might  arise  because  of  defect-Ga  complexing. 
The  results  of  these  argon  isochronal  anneals  are  shown  in 
Fig.  55.  By  800°C  all  samples  have  recovered  to  their  initial 
value  as  expected  for  these  heavily  Ga  doped  samples.  All  of 
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the  effects  seen  in  the  undoped  material  are  masked  by  the  Ga 
concentration  which  is  greater  than  lO^cm  Sample  MURR-52 
was  taken  from  the  broken  wafer.  Because  of  four  point  edge 
effects,  the  apparent  resistivity  of  this  sample  is  slightly 
higher. 

The  evaluation  of  D*  vs.  degree  of  overcompensation  has 
not  been  completed  unfortunately,  since,  at  the  close  of  this 
contract,  detectors  have  not  been  fabricated  'Tom  these  Ga 
doped  wafers  by  Rockwell  as  promised.  We  urge  that  this  set 
of  samples  be  recovered  by  AFML  for  in-house  evaluation.  We 
believe  that  the  center  sections  of  these  wafers,  at  least, 
are  suitable  for  evaluation  of  D*  vs.  Np/Ng  ratio. 

As  a final  set  of  experiments,  we  have  considered  both 
theoretically  and  experimentally  the  magnitude  of  the 
resistivity  fluctuations  as  a function  of  the  degree  of  over- 
compensation, i.e.,  the  phosphorus  to  boron  ratio.  These 
experiments  are  extremely  important  since  they  demonstrate 
conclusively  that  nuclear  transmutation  doping  provides  superior 
wafer  uniformity  for  all  phosphorus  to  boron  ratios  greater 
than  about  1.3.  For  samples  overcompensated  to  a ratio  of  2:1, 
transmutation  compensated  detector  material  possesses  resistivity 
fluctuations  which  are  about  an  order  of  magnitude  better  than 
conventionally  compensated  material. 

We  proceed  by  considering  the  fluctuations  in  con- 
ventionally compensated  material.  Typical  results4^  obtained 
at  Rockwell  for  a [P ] / [B ] ratio  of  about  1.61  are: 
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Figure  55.  Isochronal  argon  anneals  for  transmutation 
compensated  Ga  doped  Si, 
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Ap 

P 


4QQ  fi-cm 
1350  ft -cm 


29.6% 


where  we 
Ap 


will  use 

= p 

max 


the  definitions 

P • 
mm 


and 


P 


P + p . 
max  Kmin 

2 


Although  these  two  definitions  are  obviously  not  correct 

in  the  statistical  sense,  they  are  useful,  nonetheless,  in 

providing  expressions  which  can  be  easily  derived  to  provide 

considerable  insight  into  the  fluctuation  problem.  The 

conventionally  doped  experimental  data  point  is  represented 

by  the  symbol  "x"  on  Fig.  56.  It  is  not  known,  in  general, 

how  Ap/p  varies  with  p or  Np  in  conventionally  compensated 

Si;  however,  it  is  observed  from  the  NTD  literature  that 

resistivity  fluctuations  of  about  30%  have  been  observed 

47-49 

for  100  ft-cm  material  also.  We  have,  therefore,  assumed 

that  this  fluctuation  is  more  nr  less  independent  of  phosphorus 
concentration  and  have  represented  this  fact  by  the  dotted 
line  in  Fig.  56  for  the  conventional  doping  fluctuations. 

In  the  neutron  transmutation  compensated  material,  the 
only  fluctuations  in  concentration  are  those  which  were 
present  in  the  starting  material.  For  highly  zone  refined 
Si,  these  fluctuations  must  be  the  result  of  fluctuations 
in  boron  concentration.  The  phosphorus  added  by  transmutation 
is  completely  uniform,  a fact  which  will  be  demonstrated  by 
the  calculations  to  follow;  therefore,  the  impurity  fluctuation 
A[B]  is  a constant  where  [B]  represents  the  concentration  of 
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^max  - ^min 


[P]/[B] 

Figure  56.  Experimental  and  theoretical  resistivity  fluc- 
tuations as  a function  of  phosphorus  to  boron 
concentration  ratio  for  an  initial  resistivity 
fluctuation  before  compensation  of  4.721. 
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boron  in  the  starting  material.  It  is,  therefore,  obvious 
that  the  fluctuations  in  impurity  concentration  must  be 
superior  in  the  NTD  material  as  compared  to  conventional 
phosphorus  doping.  Nevertheless,  the  percentage  fluctuation 
in  resistivity  must  get  larger  near  exact  compensation  since 
very  small  changes  in  impurity  concentration  cause  rapid 
changes  in  resistivity  as  is  evident  from  Fig.  1.  Unless 
there  is  considerable  boron-phosphorus  pairing  in  convent- 
ionally doped  material,  this  fact  must  also  hold  true  in 
the  conventionally  doped  material.  Unfortunately  experiments 
to  investigate  those  fluctuations  in  conventionally  doped 
Si  have  not  been  performed. 

We  begin  an  analysis  of  the  NTD  material  with  a typical 
boron  concentration  fluctuation.  For  these  data,  we  have 
used  the  resistivity  fluctuations  in  VZ-072-#2  (Fig.  29  ) as 
representative  of  a best  case  boron  fluctuation.  From  this 
figure,  we  find  that 


. = 2.29  x 

mm 

103 

cm 

= 2.40  x 

max 

10° 

ft'Cm 

and  assuming  that  all  other  impurity  concentrations  can  be 


neglected. 

[B]  . = p . = 

L Jmm  *min 

K 

P 

253 

= 0.1054  ppb 

P • 
Kmin 

240?) 

[B]  = p = 

1 Jmax  *max 

kp 

253 

= 0.1105  pph, 

p 

pmax 

2250 

Therefore , 
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AB  = [BlmaY  - ^B]m-in  " 0.0051  ppb , 
1 Jmax  mm  rr  ’ 


[B]  = [B]  + [B]  . = 0.108  ppb, 

1 J 1 Jmax  L Jmm  " 


We  must  add  0.108  ppb  phosphorus  by  nuclear  transmutation 
for  the  average  boroi  concentration  to  be  compensated  exact- 
ly, i.e.,  [P]/[B]  = 1.  When  this  is  done,  the  highest  resis- 
tivity in  the  wafer  will  be  the  intrinsic  resistivity, 

= 214,000  Sl-cm.  Parts  of  the  wafer  will  be  n-type  and 
parts  will  be  p-type  with  resistivities  lower  than  the 
intrinsic  resistivity  due  to  impurity  fluctuations.  However, 
because  of  differences  in  the  mobility  of  electrons  and  holes, 
the  resistivity  drops  more  steeply  from  p^  on  the  p-tvpe 
side  of  intrinsic,  therefore,  the  minimum  resistivity  will 
be  in  p-tvpe  material.  We  now  have,  for  [ P ] / [ B ] = 1 

pmax  = pi  = 2*14  x 1()5 


and  we  wish  f n r a 1 r 1 1 [ 1 a t e r. 


Approximate!'7  half  of  the  car- 


rier fluctuation  will  be  on  the  p-type  side,  therefore,  the 
hole  concentration  for  this  case  can  be  calculated  from 


p = ni  + 


= 1.3  x 10ncm'3. 


0.0051 


0.0026  ppb 


The  resistivity  corresponding  to  this  concentration  is  found 
from  F.q.  (13) 


min  „ r,  2 2^ 

eMpCbni  + p ) 


and  using  y = 494cm‘”/V- sec , n^  = 1.93  x 10  cm  , and 


b = 1.6,  we  find  that 


p.  = 2.72  x 10  Q-cm 
Kmm 


Ap  = p - p • = (2.14  - 0.272)  x 10  ft-cm 

H Hmax  Hmm  v 


and 


Ap 


2.14  x 10 


5 


2.72  x 10 

rs 


= 77.91. 


2.14  x 10' 

This  calculated  point  is  represented  by  the  solid  circle  on 
Fig.  56  for  [P]/[B]  = 1.  The  second  solid  circle  at  1.3  was 
obtained  in  a similar  way. 

Since  the  number  of  intrinsic  carriers  is  negligible 
for  all  the  other  [P]/[B]  ratios  calculated,  we  will  only 
illustrate  one  of  these  calculations. 

For  [P] / [B]  = 1.63,  the  target  phosphorus  concentration 
is 

[P]  = (1.63) (0.108)  = 0.176  ppb . 

The  electron  concentration  which  is  uncompensated  is  given  by 
n = 0.176  - 0.108  = 0.068  ppb 

and 


n 

n 


87.9 


= 1292.6  Q-cm. 


n n 0.068 
To  find  the  resistivity  fluctuation,  we  need  the  electron 
concentration  fluctuation.  This  is  found  from 


n = 0.068  + -L^U- 
max  2 


= o.068  + ■Q,Q0— - = 0.0706  ppb 

2 
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and 


pmax 


Likewise , 


n 


min 


* 1245 

0.068  

0.068  - 1 .-1-  = 0.0655  ppb 


and 

pmin  0. 0655  “ 1342  n"cm* 

Therefore , 

_A£_  = 1342  - 1245  _ 7>5%  . 

p 1292.6 

This  data  point  is  represented  by  the  solid  dot  for  the 
[P]/[B]  ratio  of  1.63.  The  remainder  of  the  calculated  points 
were  obtained  similarly.  It  is  seen  that  the  fluctuations 
are  already  a factor  of  four  smaller  than  for  conventional 
doping  at  the  [ P ] / [ B ] ratio  of  1.63  and  the  improvement 
increases  for  higher  ratios.  When  the  [P]/[B]  ratio  is  3, 
thp  rpci<:ti  vi  tv  f i iirtyn t ions  become  smaller  than  the  percent 
age  error  in  the  four  point  resistivity  measurements.  We, 
therefore,  conclude  that  this  ratio  is  ideal  to  obtain  good 
uniformity  in  an  IR  detector  array.  This  ratio  is  also  close 
to  that  expected  for  a maximum  D*. 

It  should  be  noted  that  an  increase  in  resistivity 
fluctuations  in  the  starting  material  increases  the  percent- 
ages shown  in  Fig.  56  while  a decrease  would  lower  this  curve. 
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The  triangles  in  Fig.  56  represent  data  points  obtained 
experimentally  from  a set  of  similar  undoped  wafers  where 
the  [P]/[B]  ratio  was  obtained  by  nuclear  transmutation  and 
the  fluctuations  were  measured  using  four  point  resistivity 
probes.  The  data  points  at  [P]/[B]  = 4.45  and  10.0  were 
corrected  for  the  fact  that  their  initial  resistivity 
fluctuations  were  7.12%  and  5.94%  respectively  instead  of 
4.72%  assumed  in  the  model  calculation. 

The  p(x)  vs.  x data  for  the  samples  used  in  the  above 
experiment  and  taken  after  transmutation  doping  and  annealing 
are  shown  in  Fig.  57.  The  initial  resistivity  fluctuation 
and  final  [P]/[B]  ratio  is  indicated  by  each  sample  number. 

The  fluctuations  decrease  as  the  [P]/[B]  ratio  increases  in 
accordance  with  the  theory. 

We  believe  that  this  set  of  calculations  and  experiments 
demonstrates  that  NTD-Si  does  indeed  produce  a more  uniform 
resistivity  profile  and  should,  therefore,  contribute  greatly 
to  IR  detec  Lux  array  uniformity.  This  is  one  of  the  m?jr>r 
conclusions  to  be  obtained  from  this  work. 

V.  PROGRESS  AND  ACCOMPLISHMENTS 

The  object  of  this  research  program  was  to  develop  the  technique 
of  neutron  transmutation  doping  in  silicon  to  the  point  that  the 
compensation  ratio  can  be  controlled  to  a very  high  precision. 
Significant  success  in  this  task  will  allow  much  better  control 
over  the  boron  compensation  in  silicon  used  for  extrinsic  IR 
detection.  This  in  turn  should  provide  a better  control  over 
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detector  detectivity  and  temperature  dependence. 

Although  considerable  unexpected  difficulty  has  been  en- 
countered in  controlling  the  impurity  concentrations  in  silicon 
with  any  degree  of  certainty  below  0.1  ppb,  we  feel  that  real 
progress  has  been  made  in  achieving  this  goal.  It  has  also  become 
evident  that  a better  understanding  of  the  radiation  damage,  the 
oxygen  concentration  of  the  samples  and  the  interaction  between 
the  defects  and  residual  impurities  in  silicon  before  irradiation 
are  needed  in  order  to  establish  control  over  the  entire  doping 
process  in  concentrations  below  0.1  ppb. 

We  believe  that  neutron  transmutation  doping  has  removed  one 
of  the  major  difficulties  in  achieving  exact  compensation  in 
silicon,  i.e.,  the  precise  control  over  phosphorus  concentration 
which  is  improved  by  a factor  of  at  least  50:1  over  conventional 
doping  techniques.  Unfortunately,  other  residual  impurity  problems 
have  been  detected  for  the  first  time  as  a result  of  achieving  this 
control.  We  feel,  however,  that  a much  better  understanding  of 
the  limitations  involved  in  this  problem  have  resulted  from  our 
concerted  effort  on  undoped  material  where  impurity  concentrations 
can  be  determined  with  a much  higher  precision.  We  also  believe 
that  we  now  have  evidence  to  support  the  contention  that  compen- 
sation control  in  n-type  material  is  possible.  We  have  been  able 
to  demonstrate  experimentally  and  theoretically  that  NTD-Si  is 
superior  in  terms  of  wafer  uniformity  for  detector  arrays. 

The  following  is  an  abbreviated  list  of  the  significant 
accomplishments  achieved  during  this  program. 
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Construction  of  the  Bulk  Pool  Variable  Flux  Facility 
and  demonstration  of  fluence  reproducibility  to  +_  0.51. 

Demonstration  of  doping  concentration  control  to  +_  1 % 
for  100  ft-cm  material  and  +_  3%  for  resistivities  as  high  as 
100,000  fi-cm. 

No  significant  impurity  problems  detectable  by  neutron 
activation  analysis  down  to  concentrations  of  better  than  0.1 
ppb  with  the  exception  of  sodium  in  concentrations  of  1 to 
8 x 10^^atoms/cm^. 

Pre-irradiation  materials  analysis  of  p,  y,  and  t and 
activation  analysis  indicate  that  silicon  from  all  sources 
are  roughly  comparable  in  quality.  The  largest  exception  is 
in  minority  carrier  lifetime,  x,  which  is  not  as  large  as 
advertised  in  many  cases. 

Theoretical  analysis  of  the  behavior  of  n,  R^,  y,  and 
p versus  fluence  has  been  developed  to  evaluate  results  of 
neutron  doping  compensation  of  undoped  materials. 

The  theoretical  cross  section  tor  phosphorus  production 
has  been  verified  to  about  +_  1%. 

Estimates  have  been  made  on  the  number  of  displacements 
as  a function  of  fast  to  thermal  neutron  ratio.  Importance 
of  gamma  recoil  damage  calculations  have  been  verified. 

Several  differences  have  been  discovered  between  our 
annealing  results  and  those  found  in  the  literature.  These 
have  been  explained  in  terms  of  fast  to  thermal  neutron  ratio, 
and  carrier  type  immediately  after  irradiation  but  before 
annealing . 
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Have  related  in-core  1 ft-cm  radiation  and  subsequent 
annealing  results  discussed  in  literature  to  a radiation 
damage  produced  amorphous  condition  which  alters  annealing 
kinetics  drastically. 

The  activation  energy  and  frequency  factor  for  the 
electrical  activation  of  transmutation  produced  phosphorus 
during  annealing  has  been  determined  for  the  first  time. 

This  analysis  suggests  that  vacancy-phosphorus  pairs  (E- 
centers)  are  related  to  the  process  of  recovery  of  electrical 
activity. 

"Czochralski-like"  annealing  has  been  observed  in 
lightly  irradiated  silicon  (<$  'v  0.5  4>c)  which  suggests  that 
clustering  of  radiation  defects  with  residual  oxygen  is  a 
problem  which  must  be  solved  in  float  zone  doped  to  con- 
centrations below  0.1  ppb. 

The  production  of  an  acceptor  defect  in  concentrations 
between  5 x 10*^  to  5 x lO^cm  ^ by  irradiation  and  annealing 
has  beer,  demonstrated  to  be  independent  of  fluence.  This 
strongly  suggests  that  this  acceptor  is  related  to  an  in- 
trinsic residual  impurity  concentration  in  the  silicon 
before  irradiation.  We  propose  that  oxygen,  which  has  been 
decorated  with  defects  to  make  it  electrically  active,  is 
the  source  of  these  acceptors.  This  result  is  exciting  since 
it  suggests  that  oxygen  detection  in  silicon  can  easily  be 
extended  to  lower  than  present  day  detection  limits  by 
several  orders  of  magnitude. 


We  have  been  able  to  produce  a slightly  overcompensated 
sample  in  at  least  one  case  where  the  residual  uncompensated 
carrier  concentration  is  3.77  x 1010cm  3 or  about  3n^.  This 
corresponds  to  a resistivity  of  greater  than  100,000  ft-cm. 

A single  temperature  anneal  procedure  has  been  developed. 

The  radiation  effects  of  fast  neutron  fluence  have  been 
investigated. 

The  thermal  to  fast  ratio  of  position  RSAF  has  been 
measured  to  be  about  30:1. 

A series  of  15  Ga  doped  samples  have  been  compensated  to 
[P]/[B]  ratios  of  1 to  20  for  the  production  of  IR  detectors. 

Isochronal  annealing  of  Ga  doped  material  has  been 
studied. 


Extensive  theoretical  work  on  NTD  compensation  near 
intrinsic  resistivities  has  been  developed.  This  knowledge 
has  been  used  to  estimate  wafer  uniformity  as  a function  of 
phosphorus  to  boron  ratio.  It  has  been  demonstrated  con- 
clusively both  experimentally  and  theoretically  that  wafer 
uniformity  is  better  in  NTD-Si  than  conventionally  doped 
material  by  factors  of  5:1  to  30:1  over  phosphorus  to  boron 
ratios  of  2:1  to  10:1.  In  the  NTD-Si  the  resistivity 
fluctuations  for  the  ratios  decrease  from  about  +_  5%  at  a ratio 
of  2:1  to  percentages  too  small  to  measure  for  a ratio  of 
10:1.  It  is,  therefore,  apparent  that  nuclear  transmutation 
is  an  essential  technique  to  produce  IR  detector  array  uni- 
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